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1  Abstract 

Noise  equivalent  circuit  models  in  GaN  HEMTs  have  been  developed  using  both  experimental  and 
numerical  approaches.  This  work  serves  two  primary  purposes,  to  guide  device  development  by 
establishing  the  primary  causes  for  noise,  and  as  a  basis  for  models  that  will  be  used  in  transmitter 
and  receiver  design.  The  results  of  this  work  will  benefit  the  study  of  noise  mechanisms  in  other 
materials  and  devices. 

A  nonlinear  field  effect  transistor  equivalent  circuit  model  is  examined  to  identify  the  funda¬ 
mental  mechanisms  that  up-convert  baseband  1  //  noise  to  near-carrier  sideband  noise  when  the 
device  is  operated  in  the  large-signal  regime.  This  model  captures  all  physical  noise  sources  and 
nonlinearities  in  the  transistor,  and  thereby  allows  a  general  cause  and  effect  treatment.  The  noise 
sources  in  the  equivalent  circuit  model  are  determined  using  low-frequency  spectrum  analyzer  and 
microwave  noise  figure-meter  data.  Using  the  example  of  an  AlGaN/GaN  high  electron  mobil¬ 
ity  transistor,  the  developed  model  correctly  describes  both  the  measured  near-carrier  sideband 
amplitude  and  phase  noise  simultaneously. 

A  numerical  approach  to  simulate  the  intrinsic  noise  sources  within  transistors  is  described,  and 
the  impact  of  spatial  correlation  between  local  fluctuations  is  investigated.  Using  a  two-dimensional 
numerical  device  solver,  spectral  densities  for  the  gate  and  drain  noise  current  sources  and  their 
correlation  are  evaluated  using  a  Green’s  function  approach,  an  equivalent  of  Shockley’s  impedance 
field  method.  Case  studies  with  an  AlGaN/GaN  high  electron  mobility  transistor  are  supported 
by  measurement  data.  While  local  velocity  fluctuations  are  correlated,  the  assumption  of  using 
uncorrelated  diffusion  noise  source  is  found  to  be  valid  for  the  FET  geometry  studied. 
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2  Background 
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Wide  bandgap  aluminum  gallium  nitride/gallium  nitride  (AlGaN/GaN)  high  electron  mobility  tran¬ 
sistors  (HEMTs)  have  received  much  attention  for  high  frequency,  high  power  and  high  temperature 
operation,  due  to  their  high  saturated  electron  drift  velocity  (1.3 xlO7  cm/s),  large  breakdown  field 
( Ec  >  3  MV/cm),  and  good  thermal  conductivity,  particularly  with  SiC  substrates  (0  >  3.3 
W/°K-cm)  [1].  Furthermore,  the  high  speed  performance  of  AlGaN/GaN  HEMTs  can  result  in 
low  microwave  noise  characteristics,  making  these  devices  candidates  for  receivers  with  a  higher 
breakdown-related  receiver  damage  threshold.  For  example,  minimum  microwave  noise  figures  of 
0.53  dB  at  8  GHz  (100  GHz  fo  device,  Vds  =  8  V)  and  0.4  dB  at  5  GHz  (58  GHz  fo  device,  Vds  — 

1  V)  have  been  reported  for  AlGaN/GaN  HEMTs  [2,  3].  In  the  case  of  transmitter  oscillators,  low 
1//  (or  flicker)  noise  is  important  because  the  nonlinear  up-conversion  of  1//  noise  into  carrier  side¬ 
band  noise  degrades  the  system  performance.  This  low  frequency  l//-type  noise  in  AlGaN/GaN 
HEMTs  has  been  associated  with  contributions  from  the  traps  [4,  5,  6]. 

Establishing  a  noise  equivalent  circuit  model  permits  extraction  of  the  intrinsic  noise  properties 
of  AlGaN/GaN  HEMTs,  allowing  performance  prediction  and  study  of  noise  mechanisms  within 
the  device.  This  model  could  then  allow  the  prediction  of  oscillator  phase  noise,  receiver  noise, 
and  the  influence  of  physical  effects  such  as  traps  and  contacts.  Furthermore,  having  a  good  noise 
model  is  critical  for  circuit  design. 

Noise  models  can  be  broadly  classified  into  numerical  and  experimental  approaches,  both  of 
which  could  lead  to  an  equivalent  circuit  representation  with  adequate  parameterization.  Monte 
Carlo  numerical  studies  can  lead  to  port  noise  [7,  8].  A  more  tractable  Green’s  function  approach 
based  on  the  transfer  impedance  concept  of  Shockley  [9]  has  been  applied  [10].  The  experimental 
two-port  characterization  yields  three  noise  parameters,  the  minimum  noise  figure  ( Fmin ),  the  noise 
resistance  ( Rn ,  or  conductance  Gn),  and  the  optimum  source  admittance  for  minimum  noise  (Pop*  = 
Gopt  +  jBopt)  [11].  With  scattering  (S)-parameter  measurements,  these  measured  noise  parameters 
can  be  coupled  to  an  equivalent  circuit  model,  an  approach  that  is  common  for  microwave  devices 
(see,  for  example,  our  previous  work  with  noise  in  resonant  tunneling  diodes  [12]). 

Figure  1(a)  shows  the  three  physical  sources  of  intrinsic  device  noise,  where  noise  due  to  contacts 
and  parasitics  is  excluded.  The  velocity  fluctuation  noise  due  to  electron  scatter  in  the  channel 
is  independent  of  frequency  through  the  microwave  frequency  range  of  interest.  The  gate  leakage 
current  produces  shot  noise  that  can  also  be  considered  frequency  independent.  Finally,  the  capture 
and  emission  by  interface  and  surface  traps  produces  the  1//  noise.  All  these  sources  of  noise  can 
be  folded  into  the  Gaussian  noise  sources  ig  and  id  in  the  FET  equivalent  circuit  of  Fig.  1(b), 
permitting  a  frequency-dependent  representation  of  total  intrinsic  noise  (and  with  the  parasitic 
elements,  the  total  extrinsic  noise). 

There  has  been  an  evolution  of  approaches  to  develop  insight  into  the  physics  and  circuit 
elements  in  Fig.  1.  An  empirical  relation  between  the  noise  figure  and  the  circuit  parameters 
employing  fitting  parameters  was  proposed  for  GaAs  FETs  by  Fukui  [13].  More  quantitative 
equivalent  circuit  modeling  has  been  pursued  by  Statz  et  al.  [14]  and  Pospieszalski  [15],  where 
the  noise  was  represented  equivalently  as  gate  and  drain  temperature  or  resistance/conductance 
using  Nyquist  noise  sources.  In  the  case  of  the  fundamental  noise  due  to  velocity  fluctuations  of 
the  electrons  in  the  channel,  one  would  anticipate  correlation  between  these  gate  and  drain  noise 
sources.  The  extraction  of  this  correlation,  necessary  for  the  complete  frequency  dependent  model, 
has  received  little  attention.  It  has,  however,  been  predicted  numerically  [10].  While  the  correlation 
was  incorporated  into  the  convenient  extraction  model  proposed  by  Pucel  et  al.  [16],  it  was  not 
specifically  studied.  In  other  experimental  work  where  the  correlation  was  specifically  addressed, 
it  has  not  been  treated  rigorously  as  a  free  variable  [15,  17,  18].  Thus  far,  noise  equivalent  circuit 
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(a)  (b) 


Figure  1:  (a)  Fundamental  noise  mechanisms  in  AlGaN/GaN  HEMTs.  (b)  Intrinsic  noise  equivalent 
circuit  model  with  ig  and  id  representing  partially  correlated  intrinsic  noise  sources.  Inside  the 
dashed  box  is  now  a  noiseless  two-port  intrinsic  transistor,  and  vx  is  the  deterministic  drive  signal. 

models  for  GaN  HEMTs  have  received  little  attention  [3,  19,  20,  21]. 

The  features  of  GaN  HEMTs  make  them  good  candidates  as  sources,  i.e.,  for  use  in  microwave 
oscillators,  an  essential  component  in  an  RF  communication  system.  Intentionally  unstable,  oscil¬ 
lators  are  triggered  by  transients  or  noise,  and  the  amplitude  of  oscillation  increases  due  to  device 
nonlinearities  until  it  reaches  a  stable  state  with  self-limited  amplitude.  Under  large  signal  opera¬ 
tion  in  an  oscillator,  the  active  device  nonlinearities  result  in  the  up-conversion  of  the  low  frequency 
noise,  which  gives  rise  to  noise  sidebands  around  the  oscillation  frequency  in  the  spectrum,  thus  de¬ 
grading  the  spectral  purity  of  the  oscillator.  An  equivalent  circuit  model  therefore  needs  to  contain 
suitable  element  nonlinearities. 

3  Noise  in  Transistors  Under  Nonlinear  Operating  Conditions 

A  transistor  operated  under  large-signal  conditions  up-converts  the  intrinsic  1  //  noise  into  amplitude- 
modulated  and  phase-modulated  sidebands  (amplitude  and  phase  noise)  around  the  carrier,  and 
these  can  be  discriminated  in  a  measurement.  For  example,  in  GaAs  transistors,  the  measured 
amplitude  and  phase  noise  levels  have  been  found  comparable  [22].  Therefore,  a  complete  large- 
signal  noise  equivalent  circuit  should  correctly  describe  both  amplitude  and  phase  noise.  Oscil¬ 
lator  amplitude  noise  is  typically  much  smaller  than  the  oscillator  phase  noise,  due  to  the  self¬ 
amplitude-limiting  mechanism  [23].  As  a  consequence,  there  have  been  major  efforts  over  a  long 
time  period  to  understand  the  dominant  phase  noise  in  transistor-based  oscillators  (see,  for  example, 
[24,  25,  26,  27,  28]). 

A  measurement-based  nonlinear  noise  equivalent-circuit  model  allows  study  of  fundamental 
mechanisms  which  will  up-convert  the  baseband  1//  noise  to  the  near-carrier  sideband  noise.  To 
produce  this  model,  the  linear  and  nonlinear  equivalent  circuit  elements  need  to  be  determined, 
as  well  as  the  noise  sources,  i.e.,  the  baseband  1/f  noise  needs  to  be  measured.  Driving  the 
transistor  in  an  amplifier  configuration,  rather  than  as  an  oscillator,  allows  for  precise  control  over 
the  operating  conditions,  facilitating  parameterization  of  a  nonlinear  device  noise  model  [29,  30, 
31,  32].  The  resulting  nonlinear  transistor  model,  configured  as  an  amplifier,  has  been  treated 
using  harmonic  balance  simulation,  where  the  coefficients  in  a  finite  Fourier  expansion  of  voltage 
or  current  are  determined  numerically.  This  type  of  analysis  has  been  used  to  relate  a  nonlinear 
model  to  amplitude  noise  [31]  and  to  phase  noise  [29,  30].  The  impact  of  transistor  nonlinearities 
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«  on  the  baseband  noise  up-conversion  have  been  investigated  through  an  analytic  study  of  FET 
oscillators  [33,  34],  The  general  understanding  that  has  evolved  is  that  in  field  effect  transistors 
(FETs),  the  nonlinear  gate-source  capacitance  ( Cgs )  gives  rise  to  phase  noise  and  the  nonlinear 
transconductance  ( gm )  and,  to  a  lesser  degree,  the  output  conductance  (gds),  contribute  to  the 
amplitude  noise  [34].  This  view  has  led  to  simplified  nonlinear  noise  equivalent  circuits,  and  has 
also  led  to  the  use  of  only  a  gate  equivalent  noise  source  [34],  In  reality,  and  as  we  show,  this 
delineation  of  the  noise  source  contributions  is  not  possible  in  general.  Furthermore,  it  is  necessary 
to  consider  both  noise  sources  at  the  input  and  the  output. 

We  developed  a  measurement-based  large-signal  noise  equivalent-circuit  transistor  model  which 
includes  all  physical  noise  sources,  and  applied  this  model  to  an  AlGaN/GaN  high  electron  mobility 
transistor  (HEMT).  This  model  correctly  describes  both  amplitude  and  phase  noise  in  the  transistor. 
The  major  noise  sources  and  nonlinearities  are  captured  in  a  simplified  analytic  model,  thereby 
allowing  a  conceptual  theory  to  be  developed.  These  models  allow  the  fundamental  understanding 
of  how  the  near-carrier  amplitude  and  phase  noise  sidebands  are  generated  when  the  device  is 
operated  in  the  large-signal  regime.  The  AlGaN/GaN  HEMT  used  for  this  study  was  grown  on  a 
SiC  substrate  and  processed  as  described  elsewhere  [1].  The  HEMT  has  significant  1//  noise  that 
will  form  the  dominant  sideband  noise  in  a  nonlinear  amplifier  or  an  oscillator. 

The  emphasis  of  this  study  was  to  establish  the  precise  relation  between  the  intrinsic  noise 
sources  and  the  near-carrier  sideband  noise  in  terms  of  the  transistor  nonlinearities.  A  simplified 
analytical  model  for  the  in-phase  amplitude  noise  and  the  quadrature  phase  noise  provides  the  qual¬ 
itative  picture,  which  has  been  supported  with  a  more  general  nonlinear  model,  that  the  amplitude 
noise  is  dominated  by  the  nonlinear  transconductance  and  the  phase  noise  has  important  contri¬ 
bution  from  the  gate  capacitance  nonlinearity  as  well  as  the  nonlinear  transconductance.  As  the 
gate  parasitic  resistance  provides  the  vehicle  that  produces  that  phase  noise,  one  can  expect  that 
reducing  the  gate  contact  resistance  will  have  significant  impact  on  oscillator  noise  performance. 
Both  the  input  and  output  noise  sources  appear  significant  and  should  be  incorporated  into  non¬ 
linear  device  noise  models.  Although  the  nonlinear  model  was  validated  only  for  an  AlGaN/GaN 
HEMT,  it  appears  that  this  view  should  be  applicable  to  all  FETs.  The  example  large  signal 
model  used  works  effectively,  and  allowed  the  impact  of  physical  noise  sources  and  nonlinearities 
to  be  evaluated.  It  should  be  able  to  accurately  model  FETs  during  the  design  of,  for  example, 
oscillators. 

The  results  of  this  work  were  presented  at  the  2004  IEEE  International  Microwave  Sympo¬ 
sium  [32] ,  and  a  journal  paper  has  appeared  in  the  IEEE  Transactions  on  Microwave  Theory  and 
Techniques  [35].  This  journal  paper  is  included  in  the  Appendix  of  this  report. 

4  Numerical  Noise  Model 

More  complete  interrogation  of  the  device  noise  can  be  established  using  numerical  models  for  the 
noise  sources.  The  impedance  field  method  for  numerical  transistor  noise  modeling,  first  proposed 
by  Shockley  [9],  involves  the  determination  of  the  Green’s  function  (or  impedance  field),  i.e.,  the 
influence  of  a  local  fluctuation  (in  terms  of  current  or  voltage  at  one  point  in  the  channel)  at  the  gate 
and  drain  terminals,  in  the  case  of  the  field-effect  transistor  (FET).  This  method  has  enabled  the 
intrinsic  noise  sources  (gate  and  drain  noise  sources  and  their  correlation)  to  be  determined  [10,  36]. 

In  the  impedance  field  method,  the  channel  is  divided  into  many  voxels  for  numerical  simulation. 
For  long  channel  devices,  the  volume  of  each  voxel  can  be  sufficiently  large  so  that  local  velocity 
fluctuations  can  be  assumed  to  be  spatially  uncorrelated  between  voxels.  This  assumption  of 
uncorrelated  noise  sources  has  been  commonly  used  to  simplify  the  problem  and  to  calculate  the 
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intrinsic  noise  sources  [9,  10,  37].  However,  as  the  device  and  voxel  size  shrinks,  local  noise  sources 
(due  to  velocity  fluctuation)  between  voxels  can  be  correlated  because  a  significant  number  of 
electrons  can  travel  between  nearby  voxels  without  their  velocity  being  randomized  significantly. 

The  spatial  correlation  between  local  velocity  fluctuations  has  been  studied  numerically.  The 
two-point  cross-correlation  function  of  velocities  at  different  spatial  positions  can  be  related  to 
the  probability  that  electrons  travel  between  points  without  their  velocity  being  randomized  by 
scattering  events  [38] .  With  a  Monte-Carlo  simulation,  the  cross-correlation  function  of  the  velocity 
fluctuations  in  different  segments  inside  a  one-dimensional  (1-D)  GaAs  structure  was  obtained  by 
tracking  the  kinetics  of  individual  electrons  within  the  device  [39,  40].  This  demonstrated  the 
practical  importance  of  retaining  correlations  between  the  noise  sources  in  neighboring  voxels. 

We  have  used  a  two-dimensional  (2-D)  numerical  approach  to  investigate  the  effect  of  spatially 
correlated  local  noise  sources  on  the  gate  and  drain  terminal  noise  in  the  FET  geometry.  While 
Monte-Carlo  simulation  can  determine  the  numerical  values  of  the  spatially  correlated  local  noise 
sources  [39],  it  is  computationally  inefficient.  Our  proposed  approach  is  in  the  framework  of  the 
drift-diffusion  model  that  we  find  valid  for  noise  modeling  even  for  the  sub-micron  gate  length  device 
considered.  Using  this  model,  diffusion  noise  sources  and  Green’s  functions  are  determined.  We 
assume  that  the  independent  (uncorrelated)  local  diffusion  noise  sources  that  are  initially  obtained 
from  the  drift-diffusion  model  are,  in  fact,  comprised  of  contributions  from  correlated  noise  terms 
involving  the  nearby  voxels.  Then,  the  coupling  of  the  local  noise  sources  (either  correlated  or 
uncorrelated)  to  the  gate  and  drain  terminal  in  the  FET  is  simulated  using  the  Green’s  function 
approach  to  study  the  merit  of  assuming  spatially  correlated  noise  sources.  We  use  the  example  of  an 
aluminum  gallium-nitride/gallium-nitride  high  electron  mobility  transistor  (AlGaN/GaN  HEMT) 
that  was  fabricated  as  part  of  a  program  to  develop  high  power  microwave  sources  and  amplifiers  [1]. 
The  high  mobility  in  the  2-D  electron  gas  (2DEG)  not  only  results  in  high  operating  frequencies, 
but  also  low  microwave  noise,  making  AlGaN/GaN  HEMTs  suitable  for  low  noise  receivers,  hence 
the  interest  in  their  microwave  noise  properties  [41]. 

In  the  0.25  /im  AlGaN/GaN  HEMT  studied,  the  correlation  of  the  diffusion  noise  source  over 
meaningful  lengths  did  not  alter  the  terminal  noise  significantly,  leading  to  the  conclusion  that 
an  uncorrelated  model  is  sufficient.  One  would  expect  that  for  yet  smaller  device  geometries 
noise  source  correlation  would  become  more  important.  Such  smaller-scale  devices  may  require 
a  Boltzmann  transport  equation  model.  It  is  possible  that  there  is  a  regime  where  a  drift-diffusion 
noise  model  is  applicable  and  correlation  is  important,  making  the  simple  neighborhood  correlation 
concept  introduced  an  expedient  tool.  While  the  numerical  implementation  described  involved  a 
drift-diffusion  model  to  determine  both  the  noise  sources  and  the  Green’s  functions,  the  concept  of 
noise  sources  based  on  velocity  fluctuations  is  general,  as  is  the  notion  of  a  neighborhood  correlation 
scheme. 

A  paper  on  this  numerical  noise  model  was  presented  at  the  2004  IEEE  International  Microwave 
Symposium  [42] ,  and  a  journal  paper  was  submitted  to  the  IEEE  Transactions  on  Electron  Devices 
(see  the  Appendix)  [43],  which  has  now  been  revised  based  on  the  reviewers’  suggestions. 
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Abstract — A  nonlinear  field-effect  transistor  equivalent-circuit 
model  is  examined  to  identify  the  fundamental  mechanisms  that 
up-convert  baseband  1//  noise  to  near-carrier  sideband  noise 
when  the  device  is  operated  in  the  large-signal  regime.  This  model 
captures  all  physical  noise  sources  and  nonlinearities  in  the  tran¬ 
sistor,  and  thereby  allows  a  general  cause-and-effect  treatment. 
The  noise  sources  in  the  equivalent-circuit  model  are  deter¬ 
mined  using  low-frequency  spectrum  analyzer  and  microwave 
noise-figure  meter  data.  Using  the  example  of  an  AlGaN/GaN 
high  electron-mobility  transistor,  the  developed  model  correctly 
describes  both  the  measured  near-carrier  sideband  amplitude  and 
phase  noise  simultaneously. 

Index  Terms — Amplitude  noise,  gallium  nitride,  noise  measure¬ 
ment,  nonlinearities,  phase  noise,  semiconductor  device  noise. 

I.  Introduction 

TRANSMITTER  or  oscillator  noise  is  constituted  by  noise 
sidebands  modulating  a  carrier.  In  the  case  of  transistors, 
this  sideband  noise  is  primarily  attributed  to  baseband  1  // 
noise,  which  is  mixed  up  to  the  microwave  carrier  frequency 
due  to  transistor  nonlinearities.  This  near-carrier  sideband 
noise  limits  system  performance  because  it  degrades  the 
signal-to-noise  ratio  for  the  near-carrier  signals. 

A  transistor  operated  under  large-signal  conditions  up-con¬ 
verts  the  intrinsic  1//  noise  into  amplitude-  and  phase-modu¬ 
lated  sidebands  (amplitude  and  phase  noise)  around  the  carrier, 
and  these  can  be  discriminated  in  a  measurement.  For  example, 
in  GaAs  transistors,  the  measured  amplitude  and  phase-noise 
levels  have  been  found  comparable  [1],  Therefore,  a  complete 
large-signal  noise  equivalent  circuit  should  correctly  describe 
both  amplitude  and  phase  noise.  Oscillator  amplitude  noise  is 
typically  much  smaller  than  the  oscillator  phase  noise  due  to 
the  self-amplitude-limiting  mechanism  [2],  As  a  consequence, 
there  have  been  major  efforts  over  a  long  time  period  to  under¬ 
stand  the  dominant  phase  noise  in  transistor-based  oscillators 
(see,  e.g.,  [3]-[7]). 

A  measurement-based  nonlinear  noise  equivalent-circuit 
model  allows  study  of  fundamental  mechanisms  which  will 
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up-convert  the  baseband  1//  noise  to  the  near-carrier  side¬ 
band  noise.  To  produce  this  model,  the  linear  and  nonlinear 
equivalent-circuit  elements  need  to  be  determined,  as  well  as 
the  noise  sources,  i.e.,  the  baseband  l/f  noise  needs  to  be 
measured.  Driving  the  transistor  in  an  amplifier  configuration, 
rather  than  as  an  oscillator,  allows  for  precise  control  over 
the  operating  conditions,  facilitating  parameterization  of  a 
nonlinear  device  noise  model  [8]— [1 1].  The  resulting  nonlinear 
transistor  model,  configured  as  an  amplifier,  has  been  treated 
using  harmonic-balance  simulation,  where  the  coefficients  in 
a  finite  Fourier  expansion  of  voltage  or  current  are  determined 
numerically.  This  type  of  analysis  has  been  used  to  relate  a 
nonlinear  model  to  amplitude  noise  [10]  and  to  phase  noise 
[8],  [9],  The  impact  of  transistor  nonlinearities  on  the  baseband 
noise  up-conversion  have  been  investigated  through  an  analytic 
study  of  field-effect  transistor  (FET)  oscillators  [12],  [13],  The 
genera]  understanding  that  has  evolved  is  that,  in  FETs,  the 
nonlinear  gate-source  capacitance  ( Cgs )  gives  rise  to  phase 
noise  and  nonlinear  transconductance  ( gm )  and,  to  a  lesser 
degree,  the  output  conductance  (gas),  contribute  to  the  ampli¬ 
tude  noise  [13].  This  view  has  led  to  simplified  nonlinear  noise 
equivalent  circuits,  and  has  also  led  to  the  use  of  only  a  gate 
equivalent  noise  source  [13].  In  reality,  and  as  we  show,  this 
delineation  of  the  noise  source  contributions  is  not,  in  general, 
possible.  Furthermore,  it  is  necessary  to  consider  both  noise 
sources  at  the  input  and  output. 

We  present  an  example  measurement-based  large-signal 
noise  equivalent-circuit  transistor  model,  which  includes  all 
physical  noise  sources,  and  apply  this  model  to  an  AlGaN/GaN 
high  electron-mobility  transistor  (HEMT).  This  model  correctly 
describes  both  amplitude  and  phase  noise  in  the  transistor.  The 
major  noise  sources  and  nonlinearities  are  captured  in  a  simpli¬ 
fied  analytic  model,  thereby  allowing  a  conceptual  theory  to  be 
developed.  These  models  allow  the  fundamental  understanding 
of  how  the  near-carrier  amplitude  and  phase-noise  sidebands 
are  generated  when  the  device  is  operated  in  the  large-signal 
regime.  The  AlGaN/GaN  HEMT  used  for  this  study  was  grown 
on  an  SiC  substrate  and  processed  as  described  elsewhere  [14]. 
The  HEMT  has  significant  l/f  noise  that  will  form  the  domi¬ 
nant  sideband  noise  in  a  nonlinear  amplifier  or  an  oscillator. 

Section  II  describes  an  underlying  noise  theory  using  a  sim¬ 
plified  model  that  is  subsequently  related  to  the  more  general 
transistor  nonlinear  model  of  Section  III.  Noise  sideband  mea¬ 
surements  and  model  predictions  for  the  GaN  HEMT  are  pre¬ 
sented  in  Section  IV. 
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model  we  use  for  the  AlGaN/GaN  HEMT  that  can  be  applied  in 
specific  applications. 

Consider  that  vm(t)  —  v(t),  where  v(t)  is  a  deterministic 
sinusoidal  signal  given  by 

v(t )  =  Vp  cos(oj0t)  (4) 


Fig.  1.  Nonlinear  transistor  equivalent  circuit  with  capacitor  and 
transconductance  nonlinearities. 


II.  Nonlinear  Transistor  Noise  Theory 
A.  Simplified  Nonlinear  Transistor  Model 

In  the  case  of  a  nonlinear  FET,  and  considering  the  simple 
model  of  Fig.  1  (with  ac  signals  only),  the  voltage-dependent 
input  and  output  ac  currents  can  be  represented  by  a  Taylor-se- 
ries  expansion.  Assuming  weak  or  moderate  nonlinearities  and, 
hence,  terms  up  to  second  order, 

•  n\  n  dvx{t)  ,  n  dvx(t) 

^out(f)  — •  QmtVxifi)  “H  (2) 

where  vx  is  the  ac  voltage  across  the  input  capacitance,  C\  and 
gm i  are  the  linear  capacitance  and  transconductance,  and  C2 
(FA/)  and  gm 2  (A/V2)  represent  the  capacitance  and  transcon¬ 
ductance  nonlinearities,  respectively.  Equations  (1)  and  (2)  are 
coupled,  as  both  i-m  and  iout  are  a  function  of  vx .  Using  vin(t)  = 
+  vx(t)  with  (1)  leads  to 

=  RCl — +  RC2Vx(t ) — +  vx(t)  (3) 

a  nonlinear  differential  equation  in  vx.  Note  that,  in  this  simple 
model,  the  linear  series  resistor  R  is  an  important  contributor  to 
the  nonlinear  capacitance  terms.  If  R  =  0  in  (3),  v-m{t)  =  vx(t), 
and  the  capacitive  nonlinearity  will  not  have  an  impact  on  the 
output  current  *out-  The  series  resistor  R  may  correspond  to  the 
parasitic  resistance  (due  to  the  FET  gate  metal)  or  the  channel 
resistance  in  the  complete  equivalent-circuit  model. 

Note  that  a  second-order  series  expansion  for  both  iul  and 
lout  is  used  in  our  simplified  conceptual  model,  which  allows  a 
cause-and-effect  treatment  to  provide  analytical  insight.  Later, 
we  present  a  higher  order  model  for  an  AlGaN/GaN  HEMT 
equivalent  circuit,  thus  allowing  for  more  accurate  simulation, 
and  we  compare  the  harmonic-balance  simulation  results  for 
this  model  with  the  measurement  data.  In  the  actual  simulation 
and  measurement,  the  transistor  amplifier  (i.e.,  the  AlGaN/GaN 
HEMT)  is  driven  by  the  input  source  to  operate  in  the  weakly  (or 
moderately)  nonlinear  regime.  Therefore,  the  second-order  non¬ 
linearity  terms  in  the  model  dominate  those  from  higher  order 
nonlinearities,  allowing  application  of  the  simple  model  we  pre¬ 
sented  here.  Note  that,  in  the  oscillator  application,  third-order 
and  higher  order  terms  are  no  longer  negligible,  as  the  transistor 
is  typically  operated  in  a  strongly  nonlinear  regime.  While  our 
analytic  approach  assumes  a  nonlinear  model  with  terms  only 
up  to  second  order  to  investigate  the  up-conversion  mechanism, 
it  provides  a  basis  that  can  be  generalized  to  the  higher  order 


where  uj 0  is  the  angular  frequency  2irfo.  The  solution  of  (3)  for 
vx(t)  then  requires  harmonics  of  the  excitation  frequency  ujq. 
While  the  excitation  signal  is  an  even  function  in  time  (cosine), 
the  resulting  harmonics  from  (3)  for  vx(t)  should  have  both  co¬ 
sine  (even)  and  sine  (odd)  terms  due  to  the  capacitive  nonlin¬ 
earity.  Therefore,  we  assume  the  finite  Fourier  series  expansion 

N 

vx(t)  —  ao  +  ^  an  cos(nu>ot)  +  bn  sin(nwof) 

71  ~1 

N 

=  c0  +  £  Cnej™ot  +  c*ne~jnuat  (5) 

71=1 

where  an  and  bn  are  the  Fourier  coefficients  of  vx(t)  and  Co  = 
oo,  and  the  complex  Fourier  coefficient  cn  =  (1  /2)(an  —  jbn). 
The  complex  series  form  simplifies  the  calculation.  Using  N  = 
2  in  (5)  as  an  approximation  to  vx  (t),  substituting  it  into  (3),  and 
equating  terms  up  to  second  order,  we  obtain 


c0  =U 


Cl 


IL 

2 


c2  ~  -  ju0RC2 


(«' 


to; 

(7) 

(8) 


The  approximations  in  (7)  and  (8)  result  from  neglecting  terms 
under  the  assumption  that  waRC2  -C  w^RCi  <C  1,  which  is  the 
case  for  most  FETs  operating  at  microwave  frequencies.  Substi¬ 
tuting  (6)-(8)  into  (5)  gives 

V^(jJqRC2 

vx(t )  =  Vp  cos (u>0t)  +  - sin(2w0f) 


=  Vp  cos(tuo t)  -  RC2VP  cos(u>o t) 


d{Vv  cos(woZ)) 


dt 


Finally,  using  (4),  (9)  can  be  rewritten  as 


vx(t)  =  Vi n(t)  -  RC2vm(t)~^l . 

dt 


(9) 


(10) 


Applying  (10)  to  (2)  gives 

toutifi')  —  gml^init^  T  9m2‘C[n(j')  9’mlR^2'C\ni^  T7 

dt 

-2gm2RC2vfn(t)^^ 

+  W*2C«i)(^)2  (11) 

where  gm2'V\n(t)  represents  the  memoryless  nonlinearity  in  the 
transconductance,  the  third  term  on  the  right-hand  side  of  (11) 
accounts  for  the  capacitive  nonlinearity  (and,  thus,  the  memory 
effect  because  the  capacitor  is  a  charge  storage  element),  and  the 
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Ifest  two  terms  represent  coupled  transconductance  and  capaci¬ 
tance  nonlinearities.  We  note  that  the  last  two  terms  in  (1 1)  have 
not  appeared  in  prior  work  (see,  e.g.,  [13])  and  that  the  fourth 
(second  to  last)  term  in  (11)  is  necessary  to  explain  our  ampli¬ 
tude  and  phase-noise  measurements. 


B.  Baseband  Noise  Up-Conversion 

For  the  noise  random  process  n(t),  we  write  the  mean  as 
(n(i)),  where  ()  denotes  the  statistical  average,  the  autocorre¬ 
lation  as  Rnn(t,r)  (—  ( n{t)n{t  +  t))),  and  the  power  spec¬ 
tral  density  as  (|n(w)|2).  Here,  n(f)  is  either  the  voltage  across 
or  the  current  through  a  1-ft  resistor,  unless  stated  otherwise. 
Thus,  regardless  of  whether  n(t )  represents  a  voltage  or  current, 
the  normalized  power  spectral  density  is  expressed  as  (|n(o;)|2). 
We  assume  a  zero-mean  wide-sense  stationary  random  process 
n(t),  a  consequence  of  assuming  band-limited  noise,  so  that 
Rnn{r)  and  (|n(u;)|2)  form  a  Fourier  transform  pair  [15].  This 
assumption  is  valid  for  trap-related  1//  noise  sources  (even  with 
large-signal  operation)  because,  physically,  traps  with  a  long 
time  constant  cannot  follow  the  fast  carrier  signal  and  are  only 
influenced  by  the  average  (dc  bias). 

Now,  assume  that  vm(t)  —  v(t)  +  n(t),  where  n(i)  is  a 
voltage  noise  source  at  the  input  and  v(t)  is  given  by  (4).  The 
nonlinear  relation  defined  in  (1 1)  between  vx n(t)  and  is 

invariant  with  the  addition  of  the  noise  source  [13].  The  last 
term  on  the  right-hand  side  of  (11)  is  much  smaller  than  the 
other  terms  and  is,  thus,  ignored.  The  nonlinear  terms  (gm2,  C2) 
in  (11)  generate  additional  baseband  noise,  and  sideband  noise 
near  u>o  and  2aj0.  We  focus  on  the  noise  near  Wo,  no(f),  which 
is  given  by 


n0(f)  =2gm2VpCos(u}0t)n(t) 

+  gmiUaRCvVj,  sin(w0f)n(f) 

+  2gm2UoRC2Vp  sin(w0f)n2(f) 

-  gmiRC2Vpcos(uJot)^- 

-  4 gm2RC2Vp  cos(uj(d)-~~n(t). 


(12) 


Since  n0(t )  is  narrow-band,  it  can  be  represented  in  terms  of 
modulation  sidebands  as 


n0(f)  =  nj{t)  cos(w0f)  -  tcq(£)  sin(o;of)  (13) 

where  nj(f)  is  the  in-phase  (amplitude)  noise  component  and 
n<g(i)  is  the  quadrature  (phase)  noise  component,  both  with  re¬ 
spect  to  the  carrier  [16],  [17].  From  (12)  and  (13), 

ni(t)  =  2gm2Vpn(t)  -  gmlRC2Vp ^ 

-  4gm2RC2Vp^-n(t)  (14) 

nQ(0  =  -  gm\UoRC2Vpn{t)  -  2gm2WoRC2Vpn2(t). 

(15) 


Fig.  2.  Large-signal  noise  equivalent  circuit.  Inside  the  dotted  box  is  the 
intrinsic  device  with  intrinsic  noise  sources  ig  and  id-  Rg,  R,,  and  Rd  have 
extrinsic  thermal  noise,  which  are  not  shown. 


and  RC2  1,  the  last  two  terms  of  (14)  can  be  neglected.  Our 
experimental  data  indeed  indicates  that  the  contribution  from 
these  terms  to  the  overall  amplitude  noise  is  small,  and  that  nj 
can  be  reasonably  well  represented  by  the  first  term  in  (14).  In 
contrast,  our  experiment  indicates  that  both  terms  in  (15)  are 
important.  In  fact,  the  second  term  has  the  dominant  influence. 
Note  also  that  the  modulation  coefficient  (the  amount  of  up-con¬ 
version)  depends  not  only  on  the  transistor  nonlinearities  (gm2 
and  C2),  but  also  on  the  carrier  power  (through  Vp)  and  carrier 
frequency  (a>0)  in  the  case  of  the  phase  noise  t),q(£).  It  is  also 
worth  noting  in  (14)  and  (15)  that  the  capacitive  nonlinearity 
(C2)  will  not  affect  the  phase  noise  if  the  series  resistor  (a  linear 
element  in  this  model)  R  is  zero.  The  noise  dependence  on  R 
suggests  that,  if  the  series  resistance  at  the  input  can  be  low¬ 
ered  (e.g.,  the  parasitic  gate  metal  resistance),  the  near-carrier 
sideband  noise,  in  particular  the  phase  noise,  would  be  reduced. 
Also,  while  there  is  negligible  phase  noise  as  C2  — ►  0,  both  am¬ 
plitude  and  phase  noise  can  remain  as  gm2  — >  0. 

HI.  Nonlinear  Noise  Equivalent  Circuit 

As  an  example  of  a  nonlinear  model,  we  apply  the  large- 
signal  equivalent  circuit  of  Fig.  2  that  we  have  previously  found 
suitable  for  the  AlGaN/GaN  HEMT  used  in  our  noise  studies 
[18].  This  model  is  based  on  that  proposed  by  Curtice  and  Etten- 
berg  [19]  and  uses  a  nonlinear  description  of  the  drain  current, 
gate-to-source  (and  gate-to-drain)  capacitance,  and  diode  I-V 
characteristics  (Jgs  and  Iga  in  Fig.  2).  Our  prior  work  has  val¬ 
idated  this  model  based  on  comparisons  with  dc,  5-parameter, 
and  RF  power  sweep  data  [20].  The  total  (dc  +  ac)  drain-to- 
source  current  is  assumed  to  be 


Ids  —  {Ao  +  AiVi  +  A2V2  +  A3  V]3) 

•  tanh(7Vout)(l  +  AFout)  (16) 


where 


Vi  =  Vin(t  -  Tt)  (1  +  P(Vdso  -  Vout)) .  (17) 


Equation  (14)  indicates  that  the  amplitude  noise  results  from 
both  transconductance  (gm 2)  and  capacitance  (C2)  nonlinear¬ 
ities.  From  (15),  the  phase  noise  also  depends  on  both  capaci¬ 
tance  and  transconductance  nonlinearities.  If  gmiRC2  <  gm2 


Vin(i)  and  Vout(t)  are  the  gate-to-source  and  drain-to-source 
voltages  (dc  +  ac),  7  is  the  drain  current  saturation  parameter, 
A  controls  Jds  in  the  saturation  region,  rt  is  the  electron  transit 
time,  /3  describes  the  pinchoff  dependence  on  T^ut,  and  Faso  is 
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TABLE  I 

Large-Signal  Equivalent-Circuit  Parameters 
for  the  0.25  x  250  /ira2  ALGaN/GaN  HEMT 


^0 

Ai 

1  ^  | 

1  ^3 

7 

1  A 

0.1764 

0.0154 

0.7896 

P 

Tt 

Cgso 

Cgdo 

Cds 

1  Ri 

0.0128 

1.5  ps 

0.29  pF 

0.036  pF 

0.031  pF 

Vfcg 

Sex 

fc2 

cv 

fc  3 

fcA 

4.7612 

0.3780 

2.4 

0.0142 

0.6320 

Vdso 

Vg„„ 

Cpgrf 

Cpds 

Rds 

4.0  V 

-4.03  V 

0.003  pF 

0.004  pF 

0.005  pF 

337.5  n 

Rg 

Rs 

Rd 

Lg 

L. 

Ld 

6.68  n 

3.6  n 

4  n 

39  pH 

0  pH 

80.4  pH 

the  drain-bias  voltage  where  the  coefficients  A{  (i  —  0, 1, 2, 3) 
are  extracted.  The  gate-to-source  capacitance  is  described  by 


Cgs  —  Cgso 


'  1  +  /cl  tanh  (/c2(Vj in  +  Vfcg))  \ 
1  +  /cl  tanh  (fc2(VgS0  +  V/Cp))  y 


1  "h  V^ut  \ 
k  1  “h  cvV<\ so  y 


and  the  gate-to-drain  capacitance  is 


/-»  _ sy  t  /c3  tanh(/c4^ut) 

gd  gdo  x  tanh(/c4Vda<j) 


(18) 


(19) 


where  V&so  and  Vdso  are  the  bias  voltages  used  for  the  Cgso 
and  Csd0  extraction,  and  fc\,  fc2,  fc 3,  fc 4,  cv,  and  v/cg  are  fit¬ 
ting  parameters.  Table  I  lists  the  extracted  parameters  for  the 
0.25  X  250  /itn2  AlGaN/GaN  HEMT  we  studied,  obtained  by 
fitting  the  analytic  equations  (16),  (18),  and  (19)  to  the  mea¬ 
sured  dc  I-V  curves  and  the  bias-dependent  small-signal  scat¬ 
tering  parameter  data  by  means  of  a  nonlinear  curve  fit  using 
Matlab.  These  extracted  parameters  were  incorporated  into  a 
user-defined  model  in  Agilent  EEsof  ADS.1 

Equations  (2)  and  (16)  have  a  similar  mathematical  form, 
where  gm  1  and  gm 2  in  (2)  correspond  to  the  ac  components  of 
the  Ai  and  A2  terms  in  (16)  if  the  dependence  of  ids  on  Kut 
in  (16)  can  be  neglected.  Note  that  while  (1)  and  (2)  are  em¬ 
ployed  to  capture  the  basic  concept  of  transistor  nonlinearities, 
the  nonlinear  expressions  (16)  and  (18)  for  ids  and  Cgs  have 
cubic  and  hyperbolic  tangent  terms,  representing  higher  order 
nonlinearities  than  those  in  (1)  and  (2)  for  more  accurate  mod¬ 
eling.  Note  that  in  (18)  that  Cgs  is  a  function  of  not  only  Vjn, 
but  also  Kut.  Furthermore,  in  comparison  to  Fig.  1  with  the 
noise  source  (n(f))  only  at  the  input,  the  equivalent  circuit  in 
Fig.  2  includes  additional  noise  sources,  i.e.,  the  intrinsic  input 
and  output  noise  sources  ( ig ,  id)  and  the  extrinsic  device  para- 
sitics  ( Rg ,  Rs,  Rd),  which  generate  the  extrinsic  thermal  noise. 

Three  independent  noise  mechanisms  have  been  identified 
that  can  be  represented  as  equivalent  gate  and  drain  noise  cur¬ 
rent  sources,  as  shown  in  Fig.  2  [18].  The  velocity  fluctuation 
(ig,„  and  id,v)  is  due  to  the  electron  scatter  in  the  channel. 
The  1/f  noise  (ig,i/f  and  id,i/f)  results  from  the  trapping 
and  reemitting  of  electrons  due  to  surface  or  interface  states. 


'Agilent  Advanced  Design  System,  Circuit  Simulation,  Agilent  EEsof,  Agi¬ 
lent  EEsof  High-Frequency  Design  Solutions,  Palo  Alto,  CA,  2004. 
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Fig.  3.  Low-frequency  1/f  noise  measurement  setup  and  a  simplified 
equivalent  circuit  for  the  DUT  (AIGaN/GaN  HEMT)  with  the  1/f  noise 
sources  (ig,i/s  and  id,  1/ /)  at  the  input  and  output. 


The  gate  leakage  shot  noise  (*9,shot  and  id, shot)  is  associ¬ 
ated  with  the  randomness  of  the  electron  injection  into  the 
channel  over  the  gate  Schottky  barrier.  The  gate  and  drain 
noise  variance  ((|/s|2)  and  (|*d|2))  and  their  correlation  in 
terms  of  G  (iai*d)  / \/(l®sl2)  (M2))  are  incorporated  into 
the  noise  model,  and  the  values  for  each  noise  mechanism 
(velocity  fluctuation,  gate  leakage,  1/f  noise)  were  obtained  by 
using  the  measured  noise-figure  parameters  and  a  subsequent 
deembedding  routine,  a  step  taken  in  a  prior  study  [18].  The 
extrinsic  thermal  noise  (due  to  Rg,  Rs,  and  Rd)  was  modeled 
using  Nyquist’s  theorem  [(|i|2)  =  4 kT0BR,  where  k  is  the 
Boltzmann’s  constant,  To  is  the  room  temperature  (290  K),  and 
B  is  the  noise  bandwidth  (1  Hz)]  [21],  The  model  of  Fig.  2  with 
linear  elements  correctly  predicts  the  measured  frequency-de¬ 
pendent  microwave  noise  from  0.8  to  5.8  GHz  [18]. 

Under  large-signal  conditions,  the  most  critical  sideband 
noise  is  that  within  a  few  megahertz  of  the  fundamental  carrier. 
The  primary  source  of  this  noise  is  the  up-conversion  of  1  // 
noise,  which  is  the  single  dominant  noise  mechanism  in  an  FET 
at  baseband,  being  larger  than  velocity  fluctuation  noise,  gate 
leakage  current  noise,  and  thermal  noise  from  the  parasitics. 
Therefore,  accurate  determination  of  the  intrinsic  1/f  noise 
sources  below  a  few  megahertz  is  necessary  for  large-signal 
noise  modeling. 

In  order  to  determine  both  ig,i/j  and  id,i/f  and  their  cor¬ 
relation  in  the  equivalent  circuit  for  the  AlGaN/GaN  HEMT 
device-under-test  (DUT),  we  built  a  low-frequency  1/f  noise 
measurement  setup,  and  its  schematic  is  displayed  in  Fig.  3. 
In  the  measurement  frequency  range  (200  Hz-1  MHz),  most 
parasitics  in  Fig.  2  and  their  thermal  noise  contribution  can  be 
ignored,  and  the  simplified  FET  model  in  Fig.  3  is  a  reason¬ 
able  approximation.  The  noise  voltages  across  the  50-fi  gate 
termination  and  the  ac  drain  load,  i.e.,  Rds/ /Rl,  are  the  input 
and  output  noise  signals,  and  their  noise  power  per  unit  band¬ 
width  was  measured  using  a  low-frequency  spectrum  analyzer 
(HP3585A  with  a  1-Mfi  input  impedance,  which  displays  noise 
power  as  dBm/Hz).  This  data  was  converted  to  current  variance 
(|z(w)|2)(A2/Hz),  a  suitable  representation  for  the  intrinsic 
noise  sources  in  the  equivalent  circuit. 
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Frequency  (Hz) 

Fig.  4.  Measured  noise  variances  ((|i„|2),  (|«a|2),  and  (|itot|2))  for 
the  AlGaN/GaN  HEMT  from  200  Hz  to  1  MHz  with  Vda  =  3.1  V  and 
V'ga  =  —3.6  V.  The  dashed  line  is  the  estimated  (|itot|2)  using  the 
measured  (|ia|2)  and  (|f'a|2),  applied  to  (20),  neglecting  the  correlation  term 
-  The  dashed-dotted  line  is  the  measurement  system  noise 

floor. 


First,  in  order  to  obtain  the  drain  noise  current  variance 
((l*d,i//|2)).  the  shunt  noise  current  at  the  input  is  shorted 
(with  the  switch  closed  in  Fig.  3),  allowing  measurement  of 
the  drain  noise  current  variance  only,  without  the  influence  of 
the  gate  noise  current.  Second,  for  determination  of  (|*g,i//|2), 
the  gate  circuit  is  terminated  with  50  ft  (with  the  switch  open 
in  Fig.  3).  By  using  the  50-f2  termination,  a  relatively  low 
impedance  (|l/iwCgs|  »  50  fl),  most  of  *a,i//  will  flow  in 
the  50-ft  load,  and  (|*g,i//|2)  can  be  measured.  Third,  with  the 
50-17  gate  termination  (and  the  switch  open),  vx  ~  50ig!i/f. 
Thus,  the  additional  current  (50  gmig,i/f),  along  with  the  drain 
noise  current  source  *d,i//,  can  be  measured  as  a  total  output 
noise  current  variance  (|*tot|2)  given  by 

=  (50SW)2  (|ig,i//|2)  +  (|*d,i//|2) 

+  100jTO»{(i#,1//ii1//)}  (20) 

where  (|*g,i//|2),  (l*d,i//|2)>  and  (|«tot|2)  are  the  measured 
quantities,  gm  is  the  small-signal  transconductance  at  a  given 
bias  point,  and  5ft {  }  is  the  real  part.  We  used  the  gm  deter¬ 
mined  from  scattering-parameter  data  at  microwave  frequencies 
in  (20),  assuming  that  it  is  invariant  between  the  low-frequency 
and  microwave  regime. 

Fig.  4  shows  the  measured  low-frequency  (|'ig|2),  (|*d|2), 
and  (|*tot|2)  for  the  0.25  X  250  /im2  AlGaN/GaN  HEMT.  The 
resolution  bandwidth  was  set  at  100  Hz  for  the  measurement. 
The  observed  noise  is  of  a  1//“  type,  where  a  =  0.8-1 .0  under 
various  bias  conditions  over  the  entire  measurement  frequency 
range  (200  Hz-1  MHz).  As  shown  in  Fig.  4,  our  measurement 
results  for  (|jg|2),  (|*d|2),  and  (|'itot|2),  applied  to  (20),  indi¬ 
cates  that  the  correlation  term  5ft{(igii//ij  i//)}  is  negligible. 

For  the  imaginary  part  of  (*g,i //*d,i//)’  one  may  anticipate 
correlation  between  *g,i/ /  and  *d,i//  by  capacitive  coupling 
between  the  gate  and  channel  [18],  i.e.,  that  id  induces  the 
gate  noise  voltage  vg  through  transconductance,  resulting  in 


vg  =  id/ 9m  (if  perfect  correlation  can  be  assumed).  Then, 
ig  =  ju)Cssid/gm.  This  correlation  assumption  gives 


<i</>  - 

9m 


(21) 


We  have  demonstrated  experimentally  that,  in  the  case  of  ve¬ 
locity  fluctuation,  (|«d,„|2)  is  independent  of  frequency,  i.e.,  it 
is  white  noise,  and  that  (|*g,v|2)  oc  ui2  at  microwave  frequen¬ 
cies  [18].  The  bandwidth  of  the  white  noise  (|*d,„|2)  is  propor¬ 
tional  to  the  inverse  of  the  scattering  time  (typically  a  fraction 
of  a  picosecond).  The  measured  (|*g,i//|2)  in  Fig.  4  (oc  w-1)  is 
not  consistent  with  the  contribution  of  drain-induced  gate  noise 
w2C,|s(l*d,l//|2)/5m  (oc  u>)  predicted  by  (21).  This  indicates 
that  «g,i//  and  id,i/f  originate  mostly  from  independent  (thus, 
uncorrelated)  sources,  and  that  the  impact  of  any  correlation  be¬ 
tween  ig<i/f  and  id,i/f  is  not  significant.  A  similar  result  (no 
measurable  correlation  between  *g,i//  and  *d,i//)  was  found  for 
an  InP  HEMT  [22].  Therefore,  in  relation  to  the  trapping  effects 
on  1//  noise  in  AlGaN/GaN  HEMTs,  our  measurement  study 
suggests  that  the  traps  responsible  for  ig,i/f  are  distinct  from  the 
traps  responsible  for  id,i/f-  Note  that  this  is  in  contrast  to  ve¬ 
locity  fluctuation  noise,  where  our  previous  study  highlighted 
the  importance  of  correlation  [18].  It  is  likely  that  *g,i//  results 
from  traps  on  the  surface,  in  particular,  the  ungated  region  be¬ 
tween  the  gate  and  source,  and  that  id,i/f  primarily  results  from 
traps  located  near  the  channel  (e.g.,  the  AlGaN/GaN  interface). 
Prior  studies  have  shown  that  passivation  with  an  Si^N^  layer 
improves  the  large-signal  microwave  performance  [23]  and  re¬ 
sults  in  lower  1//  noise  [24],  presumably  by  suppressing  the 
surface  states. 

Note  that  stationarity  is  assumed  for  the  intrinsic  1  //  noise 
sources  that  are  incorporated  into  the  large-signal  equivalent- 
circuit  model.  This  assumption  is  valid  as  long  as  the  bias  level 
is  well  controlled,  i.e.,  fixed  in  the  amplifier  configuration,  since 
the  1/f  noise  current  variance  depends  only  on  the  dc  bias.  In 
the  oscillator  circuit,  however,  the  bias  level  is  unknown  and 
may  vary  under  strongly  nonlinear  conditions,  and  the  assump¬ 
tion  of  stationarity  has  been  revisited  to  model  the  intrinsic  1/f 
noise  sources  under  such  conditions  [25],  [26]. 

The  noise  variance  for  the  velocity  fluctuation  and  gate 
leakage  shot  noise,  obtained  under  small-signal  and,  thus,  sta¬ 
tionary  conditions,  are  used  in  the  large-signal  model  without 
modification.  The  velocity  fluctuation  and  gate  leakage  shot 
noise  are,  therefore,  assumed  to  be  band-limited,  at  least  in 
terms  of  the  model  domain.  Although  the  contribution  of  ve¬ 
locity  fluctuation  and  gate  leakage  shot  noise  to  the  near-carrier 
sideband  noise  is  negligible  under  large-signal  operation,  they 
provide  an  ultimate  limitation  On  the  microwave  noise  perfor¬ 
mance  in  the  small-signal  limit. 


IV.  Amplitude  and  Phase  Noise 

Fig.  5  shows  the  measurement  setup,  which  discriminates  be¬ 
tween  amplitude  and  phase  noise  by  adjusting  the  phase  dif¬ 
ference  between  the  local  oscillator  (LO)  and  RF  signals.  This 
measurement  procedure  employs  the  DUT  as  an  amplifier  rather 
than  as  an  oscillator  [1].  By  allowing  control  of  the  input  power 
(thus,  the  output  power),  operation  frequency,  and  bias  level,  the 
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Fig.  5.  Amplitude  and  phase-noise  measurement  setup  with  a  line  stretcher 
for  electrical  delay.  A  2-GHz  signal  from  a  low-noise  source  is  the  input  to  the 
DUT,  an  AlGaN/GaN  HEMT.  All  three  ports  of  the  power  splitter  and  mixer  are 
matched  to  50  fi. 


amplifier  configuration  for  the  DUT  allows  a  more  complete  in¬ 
vestigation  of  the  transistor  nonlinearities.  A  2-GHz  signal  from 
a  synthesizer  (HP8673G)  is  the  input  to  the  DUT,  a  0.25-mm  pe¬ 
riphery  AlGaN/GaN  HEMT.  The  output  signal  from  the  DUT 
(RF)  is  mixed  with  the  excitation  (LO),  and  the  down-converted 
spectrum  (IF)  is  measured  by  a  low-frequency  spectrum  ana¬ 
lyzer  (HP3585A,  with  a  50-0  input  impedance).  At  the  mixer, 
the  product  of  the  LO  and  RF  signals  can  be  written  as 

*if(£)  =  ((/o  +  ni(t))  cos(w0f  +  </>)-  nQ(t)  sin {wat  +  $) ) 

•ILo  cos  (tu0f)  (22) 

where  the  LO  signal  is  expressed  as  Jlo  cos(wof),  jrf  is  a 
phase-shifted  version  of  Jocos(aiof)  +  n0(f)  with  n0(t)  given 
in  (13),  and  cj>  is  the  phase  shift  between  the  RF  and  LO  signals 
that  is  achieved  by  a  line  stretcher.  The  spectrum  analyzer 
(HP3585A)  measures  the  low-pass  filtered  (frequencies  less 
than  40  MHz)  form  of  (22),  resulting  in 

,LP/^  _  ho{Io  +  ni{t))  cos  <f>  ILonQ{t)  sin  <j>  t 

®if  (t)  -  - - - - - .  (Zi) 

Equation  (23)  indicates  that  the  dc  voltage  at  the  input  of  the 
spectrum  analyzer  (monitored  by  an  HP3478A  multimeter)  is  a 
maximum  when  the  RF  and  LO  signals  are  in-phase  (cf>  —  0°) 
and  only  n/(t)  is  measured,  and  that  dc  voltage  is  zero  when 
they  are  in  quadrature  (cf)  =  90°)  and  only  nq{t)  is  measured. 
The  in-phase  and  quadrature  (I/Q)  power  spectral  densities  can 
thus  be  obtained. 

Fig.  6  shows  the  measured  and  simulated  results  for  am¬ 
plitude  and  phase  noise  (dBm/Hz)  for  the  0.25  x  250  /rm2 
AlGaN/GaN  HEMT  as  a  function  of  the  offset  frequency  from 
the  carrier.  Fig.  6  also  shows  the  measured  system  noise  floor 
(dotted  line),  summing  noise  contributions  from  every  com¬ 
ponent  used,  except  the  DUT  (the  DUT  was  replaced  with  a 
through  connection).  The  multiple  spikes  at  lower  frequencies 
result  from  the  ac  power  supply  (60  Hz)  of  the  synthesizer. 
As  shown  in  Fig.  6,  the  measured  (solid  line)  amplitude  and 
phase  noise  with  Vds  =  4  V  and  Ids  —  20  mA  both  exhibit  a 
1/ /“-type  spectrum,  where  a  ~  1  for  offset  frequencies  from 
50  Hz  to  5  kHz. 

To  validate  the  model,  the  large-signal  equivalent  circuit  of 
Fig.  2  was  incorporated  into  a  harmonic-balance  simulator, 
which  calculates  the  in-phase  ( nj )  and  quadrature-phase  noise 
(tiq).  Fig.  7  displays  the  simulation  bench  (Agilent  EEsof 
ADS  was  used),  where  a  2.4-dBm  2-GHz  signal  (the  same 
condition  as  in  the  measurement)  was  provided  as  an  input  to 


(a) 
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Fig.  6.  Measured  and  simulated:  (a)  amplitude  and  (b)  phase  noise 
for  a  0.25-mm  periphery  AlGaN/GaN  HEMT  with  Vds  =  4  V  and 

ids  =  20  mA  (Vg,  =  —4.2  V).  The  input  power  (to  the  DUT)  is  2.4  dBm 
at  2  GHz,  and  the  sideband  noise  at  the  output  load  (50  f l)  is  measured  and 
simulated  in  terms  of  power  per  normalized  bandwidth  (dBm/Hz).  The  long 
dashed  lines  are  the  simulated  case  with  <k.V/l2>  =  (l*a,i//|2>  =  0.  The 
dotted  lines  are  the  measured  system  noise  floor. 


Signal  source 


n, 

nQ 


Fig.  7.  Amplitude  and  phase-noise  simulation  bench.  The  DUT,  a  large-signal 
device  model  for  a  0.25  x  250  fim2  AlGaN/GaN  HEMT,  as  shown  in  Fig.  2,  is 
terminated  with  50  f!  at  the  input  and  output,  the  same  termination  impedance 
as  provided  in  the  actual  measurement. 


the  DUT,  and  the  I/Q  demodulator  at  the  50-0  load  produces 
two  baseband  output  signals  nj  and  r iq.  The  simulated  results 
in  Fig.  6  (square-solid  line)  show  excellent  agreement  with 
the  measured  data.  Setting  both  ig,\/f  and  in  the  model 
to  zero  produces  frequency-independent  velocity  fluctuation 
noise  and  shot  noise  shown  as  the  dashed  line  in  Fig.  6.  Under 
large-signal  conditions,  the  velocity  fluctuation  noise  is  negli¬ 
gible  at  near-carrier  sideband  frequencies.  It  is  also  shown  in 
Fig.  6  that  ig,i/f  is  a  more  important  noise  contributor  than 
id,i /f  because  turning  off  only  in  the  model  results 

in  a  greater  reduction  in  both  the  amplitude  and  phase  noise 


1320 


IEEE  TRANSACTIONS  ON  MICROWAVE  THEORY  AND  TECHNIQUES,  VOL.  53,  NO.  4,  APRIL  2005 


(a) 


(b) 

Fig.  8.  Measured  and  simulated:  (a)  amplitude  and  (b)  phase  noise.  The  input 
power  (to  the  DUT)  is  2.4  dBm  at  2  GHz.  The  squares  and  dotted  lines  are  the 
cases  with  linear  gm  and  linear  Cgs ,  respectively.  The  +  symbols  [note  that  they 
nearly  overlap  with  dotted  lines  in  (b)]  are  the  case  with  10%  of  the  estimated 
input  parasitic  resistance  R0.  The  frequency-independent  long  dashed  line  is 
the  case  for  a  linear  device. 

(dashed-dotted  lines)  than  turning  off  id,i/f  only  (diamonds) 
in  the  simulation.  However,  it  appears  that  the  contribution 
of  is  appreciable.  This  output  noise  source  appears 

as  noise  sidebands,  due  to  mixing  associated  with  either  the 
drain-source  nonlinearity  or  via  gate-drain  linear  feedback  and 
mixing  from  the  gate-source  nonlinear  elements  (it  appears 
that  the  gate-drain  nonlinearities  have  little  impact). 

After  having  established  an  accurate  nonlinear  equivalent-cir¬ 
cuit  model,  it  becomes  possible  to  evaluate  the  impact  that  each 
circuit  element  has  on  the  up-conversion.  Consider  the  impact 
of  the  nonlinearities  in  (16)  and  (18)  associated  with  depen¬ 
dencies  on  Vin  and  Vout.  Fig.  8  shows  the  simulated  amplitude 
and  phase  noise  with  the  V[n  dependency  removed  from  gm{= 
({dIds(Vin,  Kut))/9Fin)),  making  it  approximately  linear,  and 
a  linear  Cgs.  The  value  Cgs  is  set  equal  to  Cgso  by  setting  c„ 
and  /ci  equal  to  zero  in  (18),  making  it  linear.  For  gm,  ex¬ 
pressed  by  (16),  A2  and  A3  were  set  to  zero,  resulting  in  grn 
becoming  independent  of  Vin.  The  dependency  on  V^t,  which 
relates  to  the  description  of  the  nonlinear  output  conductance 
(=  ((dIds(Vin,  Vout))/dVout)),  is  relatively  weak.  As  a  conse¬ 
quence  of  these  steps,  Fig.  8(a)  shows  that  the  amplitude  noise 
is  appreciably  reduced  with  an  approximately  linear  gm,  but  not 
greatly  with  a  linear  Cgs.  On  the  other  hand,  Fig.  8(b)  shows 


that  the  phase  noise  can  be  significantly  reduced  by  linearizing 
both  gm  and  Css,  i.e.,  both  nonlinearities  are  important  con¬ 
tributors  to  the  phase  noise.  The  diamonds  in  Fig.  8,  which 
virtually  overlap  with  the  solid  lines,  are  the  case  where  Cs<\ 
was  linearized  by  setting  fc 3  equal  to  zero  in  (19).  This  simu¬ 
lation  study  indicates  that  linearization  of  Cg(j  has  little  effect 
on  the  amplitude  and  phase  noise.  With  the  removal  of  the  V\n 
dependence  from  both  gm  and  Cgs,  the  noise  in  Fig.  8  is  sub¬ 
stantially  reduced  (the  x  marks  in  Fig.  8).  Remaining  are  the 
relatively  small  contributions  due  to  the  Vout  dependency  and 
the  drain-source  and  grate-drain  nonlinearities.  Finally,  the  fre¬ 
quency-independent  long  dashed  lines  in  Fig.  8  are  the  case  for 
a  perfectly  linear  device,  where  no  1//  noise  up-conversion  oc¬ 
curs. 

Consider  again  the  simplified  analytical  model  of  Section  II, 
which  yielded  the  I/Q  noise  expressions  of  (14)  and  (15),  which 
assumed  only  an  input  equivalent  noise  source.  The  minimal  im¬ 
pact  of  the  linear  Cgs  on  the  amplitude  noise  in  Fig.  8(a)  suggests 
that  RC2  <  1  in  (14)  and  that,  indeed,  the  first  term,  which  is 
associated  only  with  the  nonlinear  transconductance,  is  domi¬ 
nant.  This  is  consistent  with  the  assumption  made  in  earlier  re¬ 
search  [13].  The  significant  impact  of  linearizing  both  gm  and 
CgS  in  the  phase  noise  of  Fig.  8(b)  suggests  that  both  terms  in 
(15)  are  significant.  Note  that  this  observation  is  in  contradic¬ 
tion  with  previous  assumptions  that  phase  noise  is  due  to  the 
nonlinear  capacitance,  and  the  contribution  from  the  nonlinear 
transconductance  can  be  neglected  [13]. 

To  evaluate  the  impact  of  reducing  the  gate  parasitic  resis¬ 
tance  Rg,  the  +  symbols  in  Fig.  8  show  the  simulated  amplitude 
and  phase-noise  performances  with  10%  of  the  extracted  gate 
parasitic  resistance.  As  indicated  in  this  figure  [and  in  (14)  and 
(15)],  the  near-carrier  sideband  noise,  in  particular,  the  phase 
noise,  can  be  reduced  by  reducing  the  gate  metal  resistance.  It 
was  found  in  our  simulation  that  lowering  the  channel  resistance 
Ri  has  little  effect  on  the  sideband  noise  (although  from  Table  I, 
both  resistances  have  similar  size).  This  is  a  consequence  of  Vin, 
the  voltage  across  Cgs  and  Ri,  modulating  Cgs  [see  (18)]  in  the 
particular  nonlinear  model  we  used.  The  point  to  note  is  the  im¬ 
pact  of  the  gate  resistance  on  the  noise  properties,  established 
through  both  the  more  general  nonlinear  model  of  Fig.  2  and  the 
simple  model  of  Fig.  1.  While  the  contribution  of  thermal  noise 
introduced  by  the  parasitics  to  the  near-carrier  sideband  noise  is 
negligible,  the  input  parasitic  resistance  has  a  significant  impact 
on  both  the  amplitude  noise  and  phase  noise  through  the  device 
nonlinearities. 

V.  Conclusion 

The  emphasis  of  this  study  was  to  establish  the  precise 
relation  between  the  intrinsic  noise  sources  and  the  near-carrier 
sideband  noise  in  terms  of  the  transistor  nonlinearities.  A  sim¬ 
plified  analytical  model  for  the  in-phase  amplitude  noise  and 
quadrature  phase  noise  provides  the  qualitative  picture,  which 
has  been  supported  with  a  more  general  nonlinear  model,  that 
the  amplitude  noise  is  dominated  by  the  nonlinear  transcon¬ 
ductance  and  the  phase  noise  has  important  contribution  from 
the  gate  capacitance  nonlinearity,  as  well  as  the  nonlinear 
transconductance.  As  the  gate  parasitic  resistance  provides  the 
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Vehicle  that  produces  that  phase  noise,  one  can  expect  that 
reducing  the  gate  contact  resistance  will  have  significant  impact 
on  oscillator  noise  performance.  Both  the  input  and  output 
noise  sources  appear  significant  and  should  be  incorporated 
into  nonlinear  device  noise  models.  Although  the  nonlinear 
model  was  validated  only  for  an  AlGaN/GaN  HEMT,  it  appears 
that  this  view  should  be  applicable  to  all  FETs.  The  example 
large-signal  model  used  works  effectively,  and  allowed  the 
impact  of  physical  noise  sources  and  nonlinearities  to  be  eval¬ 
uated.  It  should  be  able  to  accurately  model  FETs  during  the 
design  of,  for  example,  oscillators. 
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Abstract 

A  numerical  approach  to  simulate  the  intrinsic  noise  sources  within  transistors  is  described,  and  the  impact 
of  spatial  correlation  between  local  fluctuations  is  investigated.  Using  a  two-dimensional  numerical  device  solver, 
spectral  densities  for  the  gate  and  drain  noise  current  sources  and  their  correlation  are  evaluated  using  a  Green’s 
function  approach,  an  equivalent  of  Shockley’s  impedance  field  method.  Case  studies  with  an  AlGaN/GaN  high 
electron  mobility  transistor  are  supported  by  measurement  data.  While  local  velocity  fluctuations  are  correlated, 
their  impact  on  terminal  noise  is  found  identical  to  the  one  generated  by  assuming  uncorrelated  diffusion  noise 
source.  Therefore,  using  uncorrelated  local  fluctuation  to  calculate  the  intrinsic  noise  sources  is  found  to  be  valid 
even  for  the  sub-micron  gate  length  FET  studied. 
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,  I.  Introduction 

Transistor  noise  presents  a  critical  limitation,  dictating,  for  example,  receiver  sensitivity  and  transmitter  sideband 
noise.  Measurements  with  an  equivalent  circuit  model  allows  establishing  a  link  to  physical  aspects  of  the  device 
through  circuit  elements.  A  more  complete  interrogation  of  the  device  noise  can  be  established  using  numerical 
models  for  the  noise  sources.  The  impedance  field  method  for  numerical  transistor  noise  modeling,  first  proposed 
by  Shockley  [1],  involves  the  determination  of  the  Green’s  function  (or  impedance  field),  i.e.,  the  influence  of  a 
local  fluctuation  (in  terms  of  current  or  voltage  at  one  point  in  the  channel)  at  the  gate  and  drain  terminals,  in  the 
case  of  the  field-effect  transistor  (FET).  This  method  has  enabled  the  intrinsic  noise  sources  (gate  and  drain  noise 
sources  and  their  correlation)  to  be  determined  [2],  [3]. 

In  the  impedance  field  method,  the  channel  is  divided  into  many  segments  for  numerical  simulation.  For  long 
channel  devices,  the  volume  of  each  segment  can  be  sufficiently  large  so  that  local  velocity  fluctuations  can  be 
assumed  to  be  spatially  uncorrelated  between  segments.  This  assumption  of  uncorrelated  noise  sources  has  been 
commonly  used  to  simplify  the  problem  and  to  calculate  the  intrinsic  noise  sources  [1],  [2],  [4],  However,  as 
the  device  and  segment  size  shrinks,  local  noise  sources  (due  to  velocity  fluctuation)  between  segments  can  be 
correlated  because  a  significant  number  of  electrons  can  travel  between  nearby  segments  without  their  velocity 
being  randomized  significantly  [5]. 

The  spatial  correlation  between  local  velocity  fluctuations  has  been  studied  numerically.  The  two-point  cross¬ 
correlation  function  of  velocities  at  different  spatial  positions  can  be  related  to  the  probability  that  electrons  travel 
between  points  without  their  velocity  being  randomized  by  scattering  events  [6],  With  a  Monte-Carlo  simulation, 
the  cross-correlation  function  of  the  velocity  fluctuations  in  different  segments  inside  a  one-dimensional  (1-D)  GaAs 
structure  was  obtained  by  tracking  the  kinetics  of  individual  electrons  within  the  device  [7],  [8].  This  demonstrated 
the  practical  importance  of  retaining  correlations  between  the  noise  sources  in  neighboring  segments. 

In  this  paper,  a  two-dimensional  (2-D)  numerical  approach  is  presented  in  order  to  investigate  the  effect  of 
spatially  correlated  local  noise  sources  on  the  gate  and  drain  terminal  noise  in  the  FET  geometry.  While  Monte- 
Carlo  simulation  can  determine  the  numerical  values  of  the  spatially  correlated  local  noise  sources  [7],  several 
studies  have  found  that  the  assumption  of  uncorrelated  diffusion  noise  source  is  sufficient  to  determine  gate  and 
drain  noise  sources  in  the  FET  numerically  even  for  the  sub-micron  gate  length  devices  considered  [9],  [10],  The 
focus  of  our  work  is  to  investigate  how  spatial  correlation  between  local  fluctuations  affects  terminal  noise,  not 
to  determine  physical  local  noise  sources  that  are  correlated.  Our  proposed  approach  is  in  the  framework  of  the 
drift-diffusion  model  that  we  find  valid  for  noise  modeling.  Using  this  model,  diffusion  noise  sources  and  Green’s 
functions  are  determined.  We  assume  that  the  independent  (uncorrelated)  local  diffusion  noise  sources  that  are 
initially  obtained  from  the  drift-diffusion  model  are,  in  fact,  comprised  of  contributions  from  correlated  noise  terms 
involving  the  nearby  segments.  Then,  the  coupling  of  the  local  noise  sources  (either  correlated  or  uncorrelated) 
to  the  gate  and  drain  terminal  in  the  FET  is  simulated  using  the  Green’s  function  approach  to  study  the  merit  of 
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assuming  spatially  correlated  noise  sources.  We  use  the  example  of  an  aluminum  gallium-nitride/gallium-nitride 
high  electron  mobility  transistor  (AlGaN/GaN  HEMT)  that  was  fabricated  as  part  of  a  program  to  develop  high 
power  microwave  sources  and  amplifiers  [11].  The  high  mobility  in  the  2-D  electron  gas  (2DEG)  not  only  results 
in  high  operating  frequencies,  but  also  low  microwave  noise,  making  AlGaN/GaN  HEMTs  suitable  for  low  noise 
receivers,  hence  the  interest  in  their  microwave  noise  properties  [12]. 

We  develop  the  correlated  diffusion  noise  theory  in  Section  II.  The  discrete  implementation  of  the  model,  with 
a  simple  correlation  assumption  implemented,  is  described  in  Section  III.  Section  IV  applies  this  model  to  an 
AlGaN/GaN  HEMT  and  compares  the  terminal  noise  predicted  with  that  from  measurements. 

II.  Theory 

A.  Velocity  Fluctuation 

The  noise  random  processes  ni(t)  and  712(f),  where  the  subscripts  indicate  vector  component  or  position  vectors, 
are  considered  as  representing  electron  velocity,  current,  or  voltage.  We  write  the  means  as  (ni(t))  and  (712(f)), 
where  ()  denotes  the  statistical  average,  the  cross-correlation  function  as  (ni(f)ri2(f  +  r)),  and  the  cross  power 
spectral  density  as  ( nin 2),  where  *  denotes  the  complex  conjugate. 

In  the  case  of  the  FET,  the  velocity  fluctuation  is  primarily  due  to  electron  scatter  from  the  lattice  (phonons) 
and  impurities  in  the  channel.  The  random  velocity  deviation,  u(rfc,f)  at  position  vk  and  time  t,  generates  the 
zero-mean  noise  current  density  j(rfc,f)  (A/cm2)  by 

j  (r*;,  t)  =  qnku(rk,t) 

=  Alim„5ZA~9  ui(rk,t),  (1) 

Auk 

l 

where  q  is  the  absolute  value  of  electron  charge  (-1.6xl0-19  C),  nk  (cm-3)  is  the  charge  density  at  rk,  and 
the  subscript  i  denotes  the  7-th  carrier  within  the  volume  Auk  that  encloses  rk.  Using  (1),  the  two-point  (r^,  rf) 
cross-correlation  function  of  the  noise  current  density  is 

(ja(rk,  t)jp(ri,  t  +  r))  =  ^lim  ,0EE  g2  Ai/fc  Aut  (Uia(rk'  Qnjpto* t  +  T))’  (2) 

i  j 

where  the  subscripts  a  and  (3  indicate  the  projections  of  j  along  the  directions  a  and  (3.  Assuming  no  correlation 
between  different  electrons  (and  thus  their  velocities)  [6],  (2)  becomes 

( ja(rk,t)jp{rht  +  T ))  =  lim  Y]q2-^—-^—{uia{rk,t)ui0(ri,t  +  T)),  (3) 

% 

which  means  that  (ja(*k,  t)jp(*h  l  +  r))  is  nonzero  only  if  the  electrons  can  travel  from  rk  to  r  1  after  time  r  has 
passed  without  their  velocity  being  completely  randomized.  ' 

We  assume  a  zero  mean  wide-sense  stationary  random  process  u  (and  thus  j),  so  that  (uia(rk,  t)iLip(r[,  t  +  r)) 
is  independent  of  t  and  forms  a  Fourier  transform  pair  with  (uja(rfc)ti^(r/))  [13].  Taking  Fourier  transforms  of 
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both  sides  of  (3)  gives 

fehWM)  =  Jta  (4) 

i 

From  (4),  the  noise  diffusion  coefficient  Dap  can  be  defined  as  [14] 

Dap(rk,r[,u})  =  ^(uia(rk)ul0(ri)) 

1 

=  o  (uia{rk,t)uip{ri,t  +  r))e  JUJT  dr,  (5) 

"  J  ~ OO 

where  the  subscript  i  is  dropped  from  the  noise  diffusion  coefficient  since  Dap(rk,ri,u)  is  assumed  to  be  identical 
for  every  electron  within  the  volume  A uk,  observed  at  t.  For  |r|  »  ts  in  (5),  where  ts  is  the  mean-free  time  (typically 
a  fraction  of  ps)  between  collisions,  the  electron  loses  memory  of  its  initial  velocity  after  several  collisions,  i.e., 
electron  velocity  becomes  randomized,  and  thus  the  cross-correlation  function,  ( ua(t)up(t+T )),  becomes  negligible. 
The  bandwidth  of  this  random  process  is  proportional  to  the  inverse  of  the  scattering  time,  ts.  Therefore,  for 
microwave  frequencies  where  lot  • Cl  and  e_jW  ~  1, 

1  r°° 

Dap(rk,ri)  ~  -  /  (uia(rk,t)u0(ri,t  +  T))  dr,  (6) 

"  J  — oo 

which  is  frequency  independent  (i.e.,  white  noise)  in  the  microwave  regime  where  we  have  performed  numerical 
simulations.  Note  that  the  lower  frequency  (microwave)  result  for  the  noise  diffusion  coefficient  in  (6)  can  be 
demonstrated  to  be  identical  to  the  spreading  diffusion  coefficient  (Dap  =  ^(Ara(t)Arp(t)),  where  Ar  =  |rfc-rj|), 
defined  first  by  Einstein  for  Brownian  motion  [15]  and  applied  in  drift-diffusion  carrier  transport  models  [14]. 
Applying  (5)  to  (4)  gives 

i 

=  }^M2nk-^—Dap{rk,ri,u).  (7) 

Avi— >o  Avi 

Generalizing  (7)  for  vector  quantities  (u  and  j), 

(j(rfc)j*(ri))  =  lim  4g2nfc-^-D(rfc,r;,u;),  (8) 

Ai/(— >0  L\Vi 

where  (j(rfc)j*(r<))  and  D(rfc,  r^to)  are  tensors.  The  factor  of  4  (instead  of  2)  in  (8)  assumes  that  frequency 
encompasses  only  positive  values,  since  the  power  spectral  density  for  real  random  processes  is  even  in  frequency 
(i.e.,  the  Fourier  transform  assigns  half  the  power  to  negative  frequencies  and  half  to  positive  frequencies)  [13].  The 
current  density  deviation  correlation  (j(rfc)j*(r;))  in  (8)  is  proportional  to  the  spatially  correlated  diffusion  noise 
source.  Note  that  (j(rfc)j*(rj))  will  diminish  with  increasing  distance  due  to  scatter,  and  will  become  negligible  if 
the  two  points  rk  and  r /  are  separated  by  more  than  several  mean-free  paths  [4],  [7],  If  j(rfc)  and  j(r;)  are  spatially 
uncorrelated,  i.e.,  (j(rfc)j*(r/))  =  0  when  rk  ±  r (8)  can  be  simplified  to 

(j(rfc)j*(r*))  =  V«feD(rfc,r/,w)(5(rfc  -  r/),  (9) 


where  5  is  the  Dirac  delta  function. 
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B.  Non-Equilibrium  Noise 

The  current  density  correlation  of  (9)  applied  in  1-D  results  in  the  uncorrelated  current  spectral  density 

(|j(r)|2)  =  4g2n(r)D(r)<5(r).  (10) 

The  Einstein  relation  under  thermal  equilibrium  gives  ^  where  k  is  Boltzmann’s  constant,  T  is  the 

temperature,  /xo  is  the  low-field  mobility,  and  Do  is  the  spreading  diffusion  coefficient  under  equilibrium  conditions. 
Therefore,  at  equilibrium,  with  use  of  the  Einstein  relation  and  a  =  qfm,  with  o  being  the  conductivity,  (10)  becomes 

(li(r)|2)  =  4fcT(r)ff(r)tf(r), .  (11) 

where  a  —  <ro  is  the  low-field  conductivity.  The  noise  spectral  density  in  (11)  is  consistent  with  the  thermal  noise 
source  for  a  resistor  developed  by  Nyquist  [16],  where  T  —  To,  the  ambient  temperature  for  the  resistor,  and  cto  is 
the  measured  (DC)  conductivity. 

Consider  now  the  non-equilibrium  situation  where  all  electrons  are  not  at  the  lattice  temperature  and  where  there 
is  a  field-dependent  mobility  /x(D).  Equation  (11)  can  now  be  interpreted  in  terms  of  the  differential  conductance 
o(E),  which  is  equivalent  to  applying  the  Einstein  relation  with  the  interpretation  that  D  now  represents  the  small 
signal  noise  diffusion  process  [4],  This  view  is  consistent  with  the  general  velocity  correlation  interpretation  of  the 
diffusion  coefficient  in  (5).  In  applying  (1 1)  to  non-equilibrium  processes,  the  temperature  can  be  adjusted  to  produce 
the  equivalent  noise  source  (thus  T  becomes  Tn,  the  noise  temperature),  or  the  differential  noise  conductance  can 
be  used  with  T  =  To,  the  ambient  (lattice)  temperature.  Both  approaches  are  common  in  noise  equivalent  circuit 
analysis  [17],  [18].  We  take  the  latter  approach  in  a  drift-diffusion  numerical  implementation  of  noise  sources. 

C.  Green ’s  Function 

The  noise  analysis  presented  in  this  paper  is  based  on  a  Green’s  function  approach,  where,  for  the  discretized 
linear  device,  Green’s  functions  are  obtained  by  solving  Poisson’s  equation  and  the  current  continuity  equations 
using  a  2-D  numerical  device  solver  [19],  [20].  Mathematically,  Poisson’s  and  the  continuity  equations  are  coupled 
since  both  equations  include  the  charge  density.  The  linearized  equations  about  the  steady-state  solutions  (potential, 
electron  and  hole  densities)  are 

£v(r)  =  s(  r),  (12) 

where  C  is  the  linearized  differential  operator  for  the  coupled  Poisson’s  and  continuity  equations  and  u(r)  is  the 
perturbed  potential  due  to  the  small  signal  stochastic  source  term  s(r).  The  Green’s  function  for  (12)  satisfies 

CG(r,rk)  =  5(r  -  rfc),  (13) 

where  G( r,  vk)  is  the  solution  (response)  at  r  to  the  spatially  impulsive  unit  source  <5(r  —  rfc).  The  solution  of  (12) 
is  therefore 


u(r)  =  /  G(r,  rfc)s(rfc)  duk. 


(14) 


,  We  assume  that  the  source  term,  s,  results  only  from  the  noise  current  density  due  to  the  velocity  fluctuation 
of  the  electrons  in  the  channel.  Thus,  if  there  is  no  noise  current  density,  no  charge  fluctuation  occurs  (by  the 
continuity  equation),  which  translates  to  zero  perturbation  at  the  observation  point  (by  Poisson’s  equation).  We  do 
not  consider  other  intrinsic  noise  mechanisms,  for  example,  the  gate  leakage  current  shot  noise  and  the  trap-related 
1//  noise.  Since  the  current  density  appears  in  the  continuity  equation  in  the  form  of  its  divergence,  we  write  the 
source  term  s  as 


s  =  Vq.  (15) 

Applying  (15)  to  (14),  we  obtain  the  perturbed  potentials  at  the  device  terminals  as 

vp  =  j  G(rp,rfc)V  •  j(rfc)  dvk  (16) 

Vq  ~  /G(r9’r*)V  'jto)  dl>u  07) 

where  p  and  q  (or  their  location  vectors,  rp  and  rq)  denote  either  the  gate  or  drain  terminals  in  the  FET.  We  assume 
the  source  terminal  is  grounded. 

Applying  the  vector  identity  (V  ■  (a; A)  =  A  •  Va:  +  a:V  •  A)  and  the  divergence  theorem,  we  can  rewrite  (16) 
(and  likewise  (17))  as 

vp  =  ^G(rp,rk)j(rk)  ■  ds  -  J  VG(rp,  rfc)  •  j(rfc)  dvk.  (18) 

Assuming  no  surface  noise  sources,  the  surfa.ce  integral  in  (18)  vanishes,  resulting  in  [3] 

vp  =  -/  VG(rp,rfc)  -j(rfc)  dvk  (19) 

vq  =  -  J  VG(r„r,)-j(r|)  dvt.  (20) 

The  cross-correlation  voltage  spectral  density  of  the  terminal  noise  voltages  is  then 

(vpV*)  =  J  J  VG(rp,rfc)  ■  (j(rk)j*(n))  ■  VG*(r9,r, )  dvk  dut.  (21) 

Finally,  spatially  correlated  diffusion  noise  sources  in  (8)  can  be  incorporated  into  the  impedance  field  method 
formalism  in  (21)  to  calculate  the  noise  Voltage  spectral  densities  at  the  gate  and  the  drain. 


III.  Numerical  Implementation 


We  describe  a  3-D  discretization  and  then  a  2-D  implementation  with  rectangular  segments.  The  subscripts  k 
and  l  now  denote  the  /c-th  and  f-th  segments  within  the  device,  and  rk  and  r j  are  the  center  points  of  each  segment. 
We  can  thus  modify  (21)  for  the  discretized  problem  as 

(vpv*q)  =  E  E  •  VG;  AukAuh  (22) 

k  i 

where  the  correlated  diffusion  noise  source  (jy*)  from  (8)  is 
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In  the  discrete  representation,  nk,  D^,  VGpk,  and  VGqk  are  considered  constant  within  each  segment.  If  there 
is  no  spatial  correlation  of  local  velocity  fluctuations  between  segments,  i.e.,  if  D^(w)  =  0  when  k  ^  l,  we  can 
rewrite  (22)  as 

( VpV *)  =  4 q2nkWGpk  ■  D kk(u)  •  VG*fc  Auk.  (24) 

k 

In  order  to  investigate  the  impact  of  the  spatially  correlated  diffusion  noise  sources  on  the  terminal  noise,  the 
Green’s  functions  and  the  correlated  (and  uncorrelated)  diffusion  noise  sources  should  be  determined  throughout 
the  device  and  incorporated  into  (22)  and  (24)  to  compare  the  results. 

By  using  the  adjoint  network  approach  [3]  for  which  reciprocity  holds  between  the  linearized  system  (FET)  and 
its  adjoint,  the  impedance  field  G  can  be  obtained  by  applying  a  unit  current  at  the  output  terminal  p  or  q  in  the 
adjoint  network  and  finding  the  voltage  produced  at  every  internal  node  (mesh  point),  rfc,  instead  of  exciting  the 
noise  current  for  each  of  the  segments  in  the  original  system  and  observing  the  voltage  vp  or  vq  produced  at  the 
terminal.  This  results  in  significant  computational  savings  [3],  We  used  a  commercial  drift-diffusion  device  solver 
to  obtain  the  Green’s  functions  [20],  Subsequently,  using  MATLAB,  the  2-D  gradient  of  Green’s  function  was 
calculated  using  finite  differences. 

The  correlated  diffusion  noise  source  has  been  determined  by  Monte-Carlo  simulation  for  homogeneous  GaAs 
[7].  In  our  numerical  noise  study,  the  uncorrelated  diffusion  noise  source  was  determined  first,  from  which  a 
representation  for  the  correlated  diffusion  noise  source  was  deduced.  Within  the  domain  of  applicability  of  the 
drift-diffusion  model,  this  approach  allows  implementation  of  a  predetermined  spatial  correlation.  This  serves  our 
purpose  to  evaluate  the  influence  of  spatial  correlation.  To  implement  the  correct  physical  correlation  would  require, 
for  example,  a  Monte  Carlo  solution  of  the  actual  device. 

For  the  2-D  numerical  device  solution  we  employ,  it  is  assumed  that  the  Einstein  relation  holds,  and  that  the 
uncorrelated  diffusion  noise  source  is  isotropic  (i.e.,  the  tensor  D  is  isotropic).  Then,  the  diffusion  coefficient  Dk 
for  the  fc-th  segment  is  calculated  from  the  electron  mobility,  pn{E),  using  the  Einstein  relation,  assuming  a  lattice 
temperature  T.  With  isotropic  diffusivity,  (24)  can  be  rewritten  as 

(vpvI)  =  A(l2nkDk  VGpfc  •  VG*fc  Avk,  (25) 

k 

where  the  scalar  value  of  the  noise  source  Aq2nkDk  is  multiplied  with  the  dot  (inner)  product  of  the  gradients  of 
the  Green’s  functions. 

The  uncorrelated  diffusion  noise  source  (i.e.,  (jkji)  =  0  when  k  ±  l)  assumes  that  the  segment  size  is  large 
enough  (as  shown  in  Fig.  1(a))  so  that  electrons  perfectly  lose  the  memory  of  their  velocities  because  of  the  many 
scattering  within  a  segment.  However,  as  the  device  size  shrinks,  so  does  the  segment  size  (as  shown  in  Fig.  1(b)), 
and  velocity  fluctuations  in  nearby  segments  can  be  correlated.  Comparing  Figs.  1(a)  and  (b),  we  can  assume  that 
the  velocity  autocorrelation  function  for  an  electron  spatially  localized  (in  terms  of  correlation)  within  the  large 
segment  in  Fig.  1(a)  is  equal  to  the  summation  of  the  cross-correlation  functions  over  the  nearby  segments  (small 
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dashed  boxes)  in  Fig.  1(b).  The  implication  is  that  there  is  an  equivalent  diffusion  noise  source  that  represents  the 
superposition  of  a  neighborhood  of  correlated  diffusion  noise  sources,  and  in  our  work,  that  this  equivalent  noise 
source  can  be  determined  using  a  drift-diffusion  model. 

Consider  a  1-D  structure,  which  in  our  case  is  the  2DEG  at  the  AlGaN/GaN  interface.  At  the  fc-th  segment,  through 
superposition,  the  equivalent  diffusion  coefficient  can  be  written  in  general  in  terms  of  a  sum  of  cross-correlation 
terms  as  [7] 

L 

Dk  =  Yh  (26) 

i=i 

where  L  is  the  number  of  segments  (including  the  fc-th  segment)  where  spatial  correlation  persists  from  the  fc-th 
segment.  Scatter  diminishes  correlation  [7],  so  L  defines  a  neighborhood.  For  example,  if  perfect  correlation  is 
assumed  over  a  neighborhood  comprised  of  segments  within  the  mean-free  path  ( ls ),  from  (26),  Dki  —  Dk/L.  One 
could  consider  a  symmetric  neighborhood,  and  we  take  this  view  as  an  extreme  case  in  evaluating  the  impact  of 
spatial  correlation  on  the  terminal  noise.  If  the  channel  distance  variable  is  x  and  the  segment  length  is  Ax,  then 
L  —  1  +  I{2ls/Ax},  where  /{•}  is  the  integer  value.  Figure  2  illustrates  the  concept  of  how  the  spatially  correlated 
diffusion  noise  source  can  be  determined  using  this  simple  procedure.  For  example,  if  Ax  =  10  nm  and  the 
correlation  length  ls  is  20  nm,  electrons  are  assumed  to  travel  (either  left  or  right  from  the  fc-th  segment  in  Fig.  2) 
without  their  velocity  being  randomized  within  ls,  and  spatial  correlation  persists  over  5  segments  (including  the 
fc-th  segment).  Therefore,  Dki  =  Dk/ 5.  A  more  sophisticated  approach  may  involve  using  predetermined  correlation 
calculations  to  ascribe  weights  that  would  diminish  with  distance.  Also,  the  symmetry  could  be  changed  to  account 
for  the  influence  of  the  applied  electric  field,  i.e.,  electron  energy,  on  scatter  [7], 

In  our  numerical  noise  study  using  the  example  of  an  AlGaN/GaN  HEMT,  we  focused  on  the  noise  contribution 
from  the  active  channel  region  (the  2DEG).  The  2-D  numerical  solver  was  used  to  obtain  the  parameters  (i.e.,  nk, 
Dk,  Ggk,  and  Gdk)  at  all  points  within  the  device.  Considering  the  1-D  channel  only  (similar  to  the  one  shown 
in  Fig.  2),  we  determined  the  correlation  terms  Dku  assuming  (26),  and  evaluated  their  effect  in  calculating  the 
terminal  noise  using  the  impedance  field  method.  Note  that  (22)  includes  the  3-D  gradients  of  the  Green’s  functions 
and  a  tensor  Dfcj.  Because  only  the  longitudinal  segments  along  the  1-D  channel  are  considered,  we  assume  that 
Dfcj  has  the  (anisotropic)  longitudinal  diffusivity  Dki  only  for  the  correlated  diffusion  noise  source.  Therefore,  with 
longitudinal  diffusivity  ( x  component)  only,  and  applying  (23)  to  (22)  assuming  u>  =  0, 

=  E4(?2nfc^f  A^’  <2?) 

where  dGgk/dx  and  dGdi/dx  were  obtained  from  finite  differences.  As  we  will  show  from  the  numerical  simulation, 
\dG/dx\  »  \dG/dy\  for  the  Green’s  functions  involving  both  the  gate  and  the  drain.  Thus,  both  the  physical 
confinement  of  electrons  to  the  channel  and  \dG/dy\  being  small  supports  the  1-D  representation  for  noise.  The 
width  dimension  is  assumed  to  be  large  relative  to  the  (x,  y )  features,  allowing  a  2-D  numerical  simulation.  The 
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influence  of  the  correlated  diffusion  noise  on  the  calculation  of  the  intrinsic  noise  sources  can  be  investigated  by 
comparing  the  simulation  results  from  (25)  and  (27). 

IV.  Noise  Simulation 

A.  FET  Geometry 

The  numerical  device  solver  we  employed  allows  the  simulation  of  electrical  behavior  in  a  2-D  structure  assuming 
a  drift-diffusion  model  [20],  Figure  3  shows  the  2-D  device  geometry  of  an  Alo.3Gao.7N/GaN  HEMT  used  in  the 
numerical  simulation.  The  entire  region  was  divided  into  a  grid  of  mutually  orthogonal  lines,  with  solutions  to  be 
determined  at  each  of  the  mesh  points.  In  Fig.  3,  the  longitudinal  grid  spacing,  Ax,  is  10  nm.  The  vertical  grid 
spacing,  Ay,  is  1  nm,  but  is  gradually  increasing  below  the  channel.  This  device  has  a  0.25  ym  gate  length  (from 
x  =  0  to  0.25  /im  in  Fig.  3)  and  a  180  A  AlGaN  barrier  layer.  The  2-D  results  are  scaled  for  the  gate  width, 
250  ym  in  this  case,  making  this  device  consistent  with  one  for  which  we  have  experimental  noise  data  [18].  The 
gate-to-drain  spacing  (2.0  ym)  is  larger  than  the  gate-to-source  spacing  (0.5  ym)  to  alleviate  breakdown  in  the 
gate-to-drain  junction  at  higher  drain  bias  (the  device  was  fabricated  with  a  view  to  power  applications)  [11],  We 
consider  the  entire  2DEG  region  between  source  and  drain  metals  (thus  from  x  =  —0.5  to  2.25  ym  in  Fig.  3)  as 
the  intrinsic  transistor  to  be  numerically  characterized. 

While  doping  is  not  necessary  in  the  actual  device  fabrication  due  to  spontaneous  polarization  and  strain-induced 
piezoelectricity  [1 1],  the  2DEG  region  in  the  numerical  model  was  intentionally  doped  to  the  required  electron  sheet 
density.  The  doping  concentration,  the  AlGaN  barrier  thickness  (1 8  nm),  and  field-dependent  electron  mobility  model 
parameters  were  adjusted  to  match  the  measured  DC  I-V  and  AC  (S-parameter)  results. 

We  consider  mesh  points  along  the  dashed  line  (y=  19  nm)  in  Fig.  3  and  the  segments  that  enclose  each  node 
at  their  center  point.  Similarly,  in  numerical  noise  modeling  for  a  short  channel  MOSFET  using  an  uncorrelated 
diffusion  noise  source  only,  the  2-D  problem  was  simplified  to  1-D  with  the  assumption  that  the  inversion  layer  of 
the  CMOS  (active  channel)  is  very  thin  [9]. 

B.  Simulation  Results 

The  noise  simulation  results  presented  are  based  on  the  Green’s  function  approach  described  in  the  previous 
sections.  First,  the  uncorrelated  diffusion  noise  source  and  Green’s  functions  at  each  node  within  the  device  (the 
2DEG  channel  in  the  HEMT)  were  obtained  by  a  2-D  numerical  device  solution.  The  gradients  of  the  Green’s 
functions  were  calculated  by  using  MATLAB  [21],  and  the  correlated  diffusion  noise  source  was  deduced  from 
(26).  Equations  (25)  and  (27),  incorporated  into  MATLAB,  were  used  to  calculate  the  noise  voltage  spectral  densities 
at  the  gate  and  the  drain  terminals  ((|us|2),  (\vd\2),  and  (vgvd)).  The  spatially  correlated  diffusion  noise  source  Dki 
was  employed  in  (27)  to  investigate  the  impact  on  the  terminal  noise  calculation.  Finally,  we  can  transform  (|t/fl|2), 
(|i>d|2),  and  {vgv*d)  into  the  noise  current  spectral  densities,  (|?'s|2),  (\id\2),  and  (igi*d)  by  using  the  defining  relation 
between  Y-parameters  and  Z-parameters  [18],  [22], 
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.  Figure  4  shows  the  simulated  uncorrelated  diffusion  noise  source  along  the  channel,  i.e.,  the  dashed  line  (y=19  nm) 
in  Fig.  3.  The  device  is  biased  at  Vds  —  4  V  and  Ids  =  60  mA,  and  /  =  1  GHz.  As  this  bias  point  is  in  the  saturation 
region,  the  electron  charge  density  (n)  is  depleted,  which  results  in  less  diffusion  noise  (4 q2n  D)  right  under  the 
gate.  Note  that  the  ohmic  regions  (extrinsic  elements)  under  the  drain  and  source  metal  contacts  are  excluded  from 
the  noise  calculation. 

Figure  5  shows  the  Green’s  functions,  Gg(x)  and  Gd{x),  obtained  from  the  2-D  device  solution.  Note  that  Q{Gg} 
and  9{G(i},  the  imaginary  parts  of  Gg  and  Gd  obtained  at  1  GHz,  are  both  negligible  in  comparison  to  the  real 
parts  5R{Gf;}  and  51? {G^}.  Also  note  that  both  3?{GS}  and  51?{Gd}  are  varying  rapidly  in  the  channel  region  right 
under  the  gate  metal  (x  —  0  —  0.25  /jm).  Figure  6  shows  the  calculated  gradients  of  the  Green’s  functions.  The 
results  for  \dGg/dx\  and  \dGg/dy\,  and  also  \dGd/dx\  and  \dGd/dy\,  are  plotted  on  a  logarithmic  magnitude  scale 
to  more  clearly  show  the  differences.  Note  that  \dGg/dx\  »  \dGg/dy\  and  \dGd/dx\  »  \dGd/dy\,  indicating  that 
the  longitudinal  variation  of  the  Green’s  function  dominates.  As  shown  in  Figs.  5  and  6,  the  gradient  of  the  Green’s 
functions  under  the  gate-to-source  region  dominates,  and  the  gradient  of  the  drain  Green’s  function  is  an  order  of 
magnitude  larger  than  that  for  the  gate. 

We  use  the  Green’s  function  approach  to  calculate  the  gate  and  drain  terminal  noise  sources.  Figure  7  shows 
the  calculated  contribution  of  the  local  diffusion  noise  source  (uncorrelated  and  correlated)  at  the  gate  and  drain 
terminals.  The  solid  lines  are  the  cases  where  an  uncorrelated  diffusion  noise  source  is  used  for  the  gate  and  drain 
( [4q2nkDk\VGgk\ 2  and  4q2rikDk\VGdk\2,  from  (25)).  To  investigate  the  impact  of  correlation,  different  correlation 
lengths  ls  (10-50  nm)  were  chosen  for  simulation,  based  on  the  evaluation  of  mean-free  path  (average  electron 
velocity  divided  by  the  scattering  rate)  estimated  for  an  AlGaN/GaN  2DEG  [23],  [24],  The  long  dashed  lines  in 
Fig.  7  are  the  case  where  there  is  perfect  correlation  over  the  length  of  ±30  nm,  i.e.,  ls  =  30  nm,  and  Dm  =  Dyjl. 
Beyond  30  nm,  the  spatial  correlation  between  local  fluctuations  is  assumed  to  vanish.  Thus,  referring  to  (27),  the 
dashed  lines  in  Figs.  7(a)  and  (b)  were  determined  from  4g2rifc^-(^p-)  ~&jr-  and  4  q2nk^-{^£L)Yl 

Figure  7  shows  that  the  noise  contribution  of  the  correlated  diffusion  noise  source  to  the  gate  and  drain  terminal 
is  different  from  the  case  where  the  uncorrelated  diffusion  noise  source  is  used  (solid  lines).  This  difference  is 
pronounced  under  the  gate  (the  gate  metal  covers  from  x  =  0  to  0.25  fim).  This  is  because  the  diffusion  noise 
source  is  multiplied  with  the  gradient  of  the  Green’s  functions,  which  vary  appreciably  in  the  channel  region  under 
the  gate,  as  shown  in  Fig.  6.  Also  note  in  Fig.  7  that  both  the  gate  noise  and  the  drain  noise  result  mostly  from  the 
gate-to-source  region  and  under  the  gate  in  the  channel.  A  similar  trend  was  found  in  a  CMOS  device  studied  [9]. 
While  the  local  diffusion  noise  source  is  small  under  the  gate  due  to  the  charge  depletion  (see  Fig.  4),  this  noise 
source  is  multiplied  by  large  Green’s  function  gradients  and  thus  results  in  an  appreciable  noise  contribution  at  the 
gate  and  drain. 
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C.  Intrinsic  FET  Noise  Sources 


In  the  simple  FET  equivalent-circuit  model,  one  can  generally  anticipate  correlation  between  the  gate  and 
drain  noise  current  by  capacitive  ( Cgs )  coupling  [18].  For  example,  id  induces  the  gate  voltage  vg  through  the 
transconductance,  gm.  Then,  vg  —  id/ 9m  if  perfect  correlation  can  be  assumed,  and  ig  =  ju>Cgsvg  =  jtoCgsid/ gm- 
This  produces  a  noise  variance  of  (|i5|2)  =  u2C2s{\id\2)  /  g^.  We  have  demonstrated  experimentally  that,  in  the 
case  of  velocity  fluctuation,  {|*d|2)  is  independent  of  frequency  in  the  microwave  operating  range  of  the  transistor 
(white  noise),  and  hence  ( \ig\ 2)  oc  oj2,  thus  establishing  the  concept  of  capacitive  coupling  between  the  gate  and 
drain  noise  [18].  In  order  to  assess  the  effect  of  the  correlation  between  ig  and  id,  the  correlation  coefficient,  C,  is 
defined  in  normalized  form  as  [18],  [25] 


C  = 


(jgjd) 

V(M2)(M2)' 


(28) 


With  perfect  correlation  assumed  in  the  ideal  noise  model,  C  —  jl. 

Figures  8  and  9  show  the  simulated  (|i9|2),  (|id|2),  and  the  correlation  coefficient  C  (magnitude  and  phase), 
transformed  from  (|ufl|2),  (|v<f|2),  and  (vgvd)  that  are  obtained  by  (25)  and  (27).  We  compare  these  results  with 
those  extracted  from  a  noise  measurement  of  the  0.25x250  gm2  Alo.3Gao.7N/GaN  HEMT  [18].  Both  measurements 
and  simulations  result  from  velocity  fluctuation  only,  with  the  influence  of  shot,  1//,  and  extrinsic  thermal  noise 
deembedded  [18].  The  drain  noise  simulation  almost  precisely  matches  the  experimental  result.  The  discrepancy 
in  the  gate  noise  result  is  due  to  the  Schottky  contact  attributes  not  being  correctly  implemented  in  the  numerical 
device  solver  [20],  We  also  noticed  that  the  gate  Green’s  function  was  sensitive  to  the  AlGaN  barrier  layer  thickness 
(nominally  20  nm),  which  might  be  different  from  that  in  the  actual  device  fabricated.  The  correlation  coefficient 
trend  is  also  captured  in  the  simulation  result  of  Fig.  9.  Using  the  diffusion  noise  source  without  correlation  and 
assuming  perfect  correlation  over  up  to  30  nm  produce  virtually  the  same  noise  results.  This  occurs  because  the 
noise  sources,  multiplied  with  the  gradients  of  the  Green’s  functions,  are  integrated  over  space,  which  tends  to 
wash  out  the  features  evident  in  Fig.  7.  Note  in  Figs.  8  and  9  that  assuming  perfect  spatial  correlation  over  ls  = 
50  nm  for  the  diffusion  noise  source  generates  slightly  different  results  for  the  terminal  noise  (especially  (|i<i|2)  and 
|C|).  However,  it  is  our  understanding  that  50  nm  is  several  mean-free  paths,  and  thus  electron  velocities  should 
be  randomized  after  traveling  such  distance. 

With  the  validated  drift-diffusion  noise  model,  Fig.  10  shows  the  simulated  (|i<i|2)  for  different  bias  conditions. 
Note  that  (|id|2)  increases  with  Id.s,  but  does  not  change  appreciably  with  Vds  in  the  saturation  regime,  as  shown  in 
Fig.  10.  This  observation  is  consistent  with  those  from  our  measurements  [18].  The  current-dependent  noise  arises 
from  non-equilibrium  transport  in  the  channel,  with  both  noise  sources  and  Green’s  functions  being  a  function  of 
bias. 

To  first  order,  high  channel  mobility  and  high  unity  current  gain  frequency  fx  produce  low  velocity  fluctuation 
noise.  The  2DEG  mobility  for  an  AlGaN/GaN  HEMT  has  been  measured  to  be  1 500  cm2/V-sec  at  300  K  and  7900 
cm2/V-sec  at  80  K  [12].  The  device  we  made  measurements  on  had  fx  =  40  GHz.  The  noise  in  this  relatively 
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high-speed  device  is  modeled  quite  nicely  using  a  drift-diffusion  model. 

V.  Conclusion 

Numerical  noise  modeling  allows  interrogation  of  intrinsic  properties  that  contribute  to  measured  terminal  noise. 
While  Monte  Carlo  techniques  provide  a  more  rigorous  physical  model,  the  simple  drift-diffusion  model  used 
predicts  the  velocity  fluctuation  noise  in  the  0.25  pm  AlGaN/GaN  HEMT  studied.  The  correlation  of  the  diffusion 
noise  source  over  meaningful  lengths  did  not  alter  the  terminal  noise  significantly,  leading  to  the  conclusion  that  an 
uncorrelated  model  is  sufficient.  One  would  expect  that  for  yet  smaller  device  geometries  noise  source  correlation 
would  become  more  important.  Such  smaller-scale  devices  may  require  a  Boltzmann  transport  equation  model.  It  is 
possible  that  there  is  a  regime  where  a  drift-diffusion  noise  model  is  applicable  and  correlation  is  important,  making 
the  simple  neighborhood  correlation  concept  introduced  an  expedient  tool.  While  the  numerical  implementation 
described  involved  a  drift-diffusion  model  to  determine  both  the  noise  sources  and  the  Green’s  functions,  the 
concept  of  noise  sources  based  on  velocity  fluctuations  is  general,  as  is  the  notion  of  a  neighborhood  correlation 
scheme.  J 
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Figure  Captions 

Fig.  1.  (a)  A  large  segment  (dashed  box)  for  which  the  uncorrelated  diffusion  noise  source  may  be  an  adequate 
representation,  (b)  Smaller  segments  where  the  electron  velocity  fluctuations  between  segments  (smaller  square 
dashed  boxes)  are  spatially  correlated. 

Fig.  2.  Discretized  1-D  structure.  The  diffusion  noise  source  for  the  fc-th  segment  has  noise  contributions  over 
neighboring  segments,  with  their  cross-correlation  terms  (denoted  as  D^i)  dependent  on  how  long  the  correlation 
length  is. 

Fig.  3.  AlGaN/GaN  HEMT  device  geometry  used  in  the  numerical  simulation.  The  longitudinal  (along  with  the 
channel)  and  vertical  grid  spacings  are  represented  as  Ax  and  Ay,  respectively.  The  gate  length  is  0.25  p,m. 

Fig.  4.  Simulated  local  (uncorrelated)  diffusion  noise  source  along  the  2DEG  channel  for  the  AlGaN/GaN  HEMT.  >> 

The  device  is  biased  at  Vds  =  4  V,  Ids  =  60  mA.  /  =  1  GHz.  The  AlGaN/GaN  HEMT  has  0.25  /im  gate  length 
(from  x  =  0  to  0.25  p. m). 

Fig.  5.  Simulated  Green’s  function  (impedance  field)  for  (a)  the  gate  and  (b)  the  drain.  The  device  is  biased  at  Vds 
=  4  V,  Ids  =  60  mA,  and  /  =  1  GHz.  The  complex  Green’s  function  has  units  f2,  and  5ft {  }  and  5{  }  are  the  real 
and  imaginary  parts. 

Fig.  6.  Calculated  magnitude  of  the  gradient  of  the  complex  Green’s  function  (impedance  field)  for  (a)  the  gate 
and  (b)  the  drain.  The  device  is  biased  at  Vis  =  4  V,  Ids  =  60  mA,  and  /  =  1  GHz. 

Fig.  7.  The  contribution  of  the  local  diffusion  noise  source  (noise  source  multiplied  with  the  gradient  of  the  Green’s 
function  and  its  conjugate)  to  (a)  the  gate  and  (b)  the  drain  noise.  The  device  is  biased  at  Vds  =  4  V,  Ids  =  60  mA, 
and  /  =  1  GHz. 

Fig.  8.  Simulated  (a)  ( \ig\ 2)  and  (b)  (\id\2)-  The  device  is  biased  at  Vds  =  4  V,  Ids  =  60  mA.  The  symbol  (x)  is 
the  measurement  result  (velocity  fluctuation  only,  after  deembedding  shot  and  1//  noise)  from  the  0.25x250  /im2 
AlGaN/GaN  HEMT. 

Fig.  9.  (a)  The  magnitude  and  (b)  phase  of  the  simulated  correlation  coefficient.  The  device  is  biased  at  Vds  =  4  V, 

Ids  =  60  mA.  The  symbol  (x)  is  the  measurement  result  (velocity  fluctuation  only,  after  deembedding  shot  and  1// 
noise)  from  the  0.25x250  /im2  AlGaN/GaN  HEMT. 

Fig.  10.  Simulated  (|Zd|2)  (velocity  fluctuation  only).  The  device  is  biased  at  Ids  =  30,  60,  90  mA,  and  for  Ids  = 

60  mA,  Vds  =  4,  6,  8  V. 
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